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SECTION  I 
INTRODUCTION 


The  purpose  of  this  manual  Is  to  provide  the  EXEMPT  user  with  an 
understanding  of  the  capability  of  the  EXEMPT  code  and  the  Information 
necessary  to  correctly  Implement  the  program. 

The  remainder  of  this  section  discusses  the  EXEMPT  code  at  a 
general  level.  Subsequent  sections  discuss  the  EXEMPT  Command  Language  ' 

(ECL),  the  external  file  structures,  the  Internal  data  structures,  the 
user  interface  requirements,  and  the  computer  requirements. 

1.  FUNCTION  OF  EXEMPT 

\ The  function  which  the  EXEMPT  code  is  designed  to  perform  is  that 

of  a data  base  interface  processor.  By  manipulating  the  data  provided 
by  the  user  and  generated  Internally,  a number  of  mathematical  operations 
may  be  performed  on  the  data.  The  operations  are  commanded  by  the  user 
and  are  performed  in  the  sequence  specified.  The  result  of  the  operations 
will  be  the  analysis  desired  on  the  system  to  which  the  data  pertains. 

Therefore,  EXEMPT  stores  and  retrieves  data,  provides  data  arrays  to 
user  supplied  programs  and  subroutines,  performs  algebraic  operations  on 
data,  and  will  perform  EXEMPT  library  functions  on  data. 

2.  DATA  SET  DESCRIPTION 

The  basic  element  of  EXEMPT  is  the  data  set.  This  is  a collection 
of  numbers  which  have  some  physical  or  logical  reason  to  be  included  in 
a set.  In  order  to  assure  the  integrity  of  the  data  and  the  sanity  of 
an  operation,  a certain  amount  of  information  about  the  data  set  is 
required  and  will  be  referred  to  as  data  set  attributes.  The  attributes 
of  a data  set  are  used  to  identify  the  type,  source,  length,  and  location 
of  each  data  set  referenced.  The  user  may  reference  these  data  sets 
using  symbolic  names  (limited  to  ten  characters).  When  a reference  is 
made  to  perform  an  operation  on  a data  set,  the  attributes  will  be 
recovered  to  determine  the  correctness  of  the  operation,  auxiliary  data 
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sets  required,  and  the  appropriate  input/output  (I/O)  action  to  recover 
the  data  and  store  current  data.  EXEMPT  will  distinguish  between  three 
basic  types  of  data  sets.  These  are  frequency  domain,  time  domain,  and 
raw  or  untyped  data.  Associated  with  frequency  domain  data  are  two 
subsets,  the  dependent  or  resoonse  data  and  the  independent  or  frequency 
data.  Similarly,  a time  domain  data  set  will  have  dependent  or  response 
data  and  independent  or  time  data.  Raw  data  constitute  a user  provided 
array  of  data  which  is  not  a frequency  or  time  domain  data  set.  Raw 
data  may  be  used  in  direct  algebraic  manipulations  as  an  arithmetic 
operand  but  may  not  have  frequency  or  time  domain  operations  performed 
upon  it.  Data  sets  provided  by  the  user  from  some  external  source  must 
have  the  type  specified.  Data  sets  generated  internally  by  operations 
on  existing  data  sets  will  have  the  same  type  attributes  as  the  existing 
or  parent  data  sets  as  modified  by  the  operation.  For  instance,  if  a 
data  set  is  generated  by  performing  a Fourier  transform  upon  a time 
domain  data  set,  the  resultant  data  set  will  have  the  attributes  of  a 
frequency  domain  data  set.  However,  in  the  direct  point  by  point  mul- 
tiplication of  two  like  data  sets,  or  a raw  and  typed  data  set,  the 
result  will  be  a data  set  of  the  same  type  as  involved  in  the  operands. 
In  order  to  avoid  unnecessary  storage,  the  independent  array  is  not 
regenerated  for  the  resultant  array  unless  the  operation  alters  the 
array  contents.  Instead,  the  resultant  dependent  data  array  is  linked 
to  the  independent  data  array  of  its  parents.  This  type  of  implementa- 
tion results  in  EXEMPT  being  an  attribute  linked/ set  theoretic  data  base 
manipulator.  As  such,  it  is  easily  expanded  to  include  other  operations 
by  providing  the  attributes  required  for  the  operations  and  the  software 
to  implement  the  operation. 


SECTION  II 

EXEMPT  COMMAND  LANGUAGE 


EXEMPT  is  a highly  modularized  code  with  distinct  classes  of  functions 
performed  in  submodules.  The  top  level  structure  of  EXEMPT  is  shown  in 
Figures  1 and  2.  As  illustrated,  all  codes  required  for  data  inter- 
facing with  system  devices  (tape  files,  disk  files,  and  punch  cards)  are 
located  in  the  main  segment.  The  first  level  of  modularization  is 
represented  by  the  Language  Input  Processor  (LIP)  and  the  Command 
Language  Processor  (CLP) . These  two  modules  execute  independently  of  each 
other.  The  LIP  is  called  to  read  the  users  ECL  and  generate  CLP  decodable 
tasks  in  the  task  table.  After  encountering  an  END  or  RUN  command, 
control  is  returned  to  EXEMPT  which  subsequently  causes  the  CLP  to  be 
loaded.  The  structure  of  the  CLP  is  shown  in  Figure  2.  There  are  three 
classes  of  operations  illustrated.  The  logical  operations  are  I/O 
oriented.  They  simply  format  and  transfer  data  to  the  appropriate 
device.  The  SLAM  operations  generate  and  modify  data  as  required  by  the 
particular  submodel  being  used.  The  arithmetic  operations  are  the 
primary  operations  within  EXEMPT.  These  operations  are  employed  in 
performing  the  analysis  desired. 

The  following  discussion  of  the  ECL  is  presented  by  module  and 
submodule.  Each  command  is  presented  on  a separate  page  to  facilitate 
use  of  the  manual  and  the  continued  development  of  EXEMPT.  All  commands 
are  input  in  a free  field  format  with  the  command  name  specified  as 
required  by  the  particular  command. 

The  commands  of  the  language  are  listed  in  Tables  1 and  2.  Each 
command  is  separated  by  at  least  one  blank  from  its  parameter  field  and 
may  start  at  any  column  in  the  scan  field.  The  parameters  for  each 
command  are  separated  by  commas  or  blanks.  There  may  be  multiple  commands 
on  a single  card.  Multiple  commands  must  be  separated  by  a double  slash 
(//). 

Continuation  is  implied  by  a termination  with  a comma  or  an  arith- 
metic operator  (+,  -,  /,  *,  ■) . Comments  are  allowed  by  placing  a S in 
any  column  and  the  comment  following  the  $.  There  is  no  continuation 
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Figure  1.  EXEMPT  Program  Structure 


COMMAND 


Figure  2.  EXEftPT  CLP  Structure 


Table  1 

EXEMPT  COMMAND  REWORDS 


ABS 

MIN 

ALOG 

MOD 

ALTFIL 

PLOT 

ATAN 

POLY 

ATAN2 

PRTPLT 

CALL 

PUNCH 

CLEAN 

PURGE 

CMP  LX 

READ 

CONJG 

READC 

COS 

READF 

DELETE 

REWIND 

DROP 

RUN 

EDIT 

SAVE 

END 

SAVRAN 

EOF 

SCAN 

EXP 

SHIFT 

FILE 

SIGN 

IORMAT 

SIN 

FT 

SKIPB 

GETRAN 

SKIPF 

INTERP 

SQRT 

INVFT 

TAN 

LABEL 

TITLE 

LIST 

TRDEXP 

LOOP 

TRDSIN 

MAX 

TYPE 

MERGE 

WRITE 
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Table  2 

SLAM  COMMAND  KEYWORDS 


AS1918 

AVBAY 

AVGSKT 

COAX 

COCKPT 

CO’^'JIT 

COUPLER 

DBCAB 

DBCAB1 

ENGINE 

HFANT 

INDCOU 

MARBEA 


PITOT 

PNLJNT 

RAD ALT 

RADAR 

RVBEA 

SFIELD 

SHIELD 

SNGLIN 

WEABAY 

WEBAOP 

WROOT 

WWELL 

WWGND 


The  following  character  set  is  recognized  by  EXEMPT. 


All  alpha  characters 

All  number  characters 

Plus  sign 

Minus  sign 

Multiply  sign 

Divide  sign 

Left  parenthesis 

Right  Parenthesis 

Dollar  sign 

Equal  sign 

Blank 

Comma 

Period 


A number  may  be  an  integer,  a real  number,  or  a number  in  E format, 
and  may  contain  a unary  operator. 

Defaults  for  missing  parameters  are  described  where  applicable; 
otherwise,  the  parameter  must  be  present. 

The  symbols  on  the  commands  are: 


l ] 

Multiple  choice  of  required  field 

Optional  fields 

Name 

User  defined  symbolic  name 

N 

Integer  number 

V, Value 

Real  number 
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1.  LIP  COMMANDS 

In  order  to  provide  the  user  additional  flexibility  during  LIP 
execution,  four  commands  are  recognizable  by  the  LIP.  These  commands 
must  be  the  only  commands  on  the  card;  however,  they  may  appear  anywhere 
on  the  card.  They  are  used  by  the  prescanning  program  which  collects 
all  of  the  user's  input  (ECL  and  coded  data)  from  the  card  input  files 
and  user  specified  alternate  input  file  and  rewrites  the  data  onto  a 
file  to  be  read  by  the  LIP.  These  commands  will  not  be  transferred  to 
the  CLP.  The  normal  input  file  is  FORTRAN  logical  unit  5 and  the  alter- 
nate file  available  is  FORTRAN  logical  unit  4.  If  used,  the  user  must 
attach  this  file  as  TAPE4  and  it  must  be  a coded  card  image  file. 

EXEMPT  reads  all  coded  data  from  these  files  and  writes  the  data  on 
FORTRAN  logical  unit  3.  Therefore,  the  user  should  avoid  using  a local 
file  name  of  TAPE3  unless  he  desires  to  save  the  EXEMPT  input  data  after 
initial  processing. 


; 
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a.  Scan  Limit  Directive 


SCAN 


SCAN  [Nl]  ( ,N2] 

This  directive  establishes  the  columns  of  the  input  cards  to 
be  searched  for  ECL  directives.  Nl  specifies  the  first  column  and  if 
omitted  is  assumed  to  be  1.  N2  specifies  the  last  column  and  is  not 
defaulted.  If  a SCAN  command  is  not  provided,  card  columns  1-80  are 
used.  If  Nl  is  not  provided,  a comma  must  precede  N2.  All  data  preceding 
column  Nl  and  following  column  N2  are  ignored.  If  an  END  or  RUN  command 
is  encountered,  the  scan  limits  will  be  reset  to  columns  1-80  in  order 
to  read  data.  Thus,  a SCAN  directive  is  required  after  each  RUN  command 
to  reestablish  the  desired  scan  limits. 

Examples: 

Establish  the  scan  limit  of  columns  1 to  50. 

SCAN  ,50 

Establish  the  scan  limit  of  columns  10  to  60  with  RUN  directive. 

SCAN  10,  60 


ECL 


RUN 


n 


: 


CARD  DATA 


END  OF  DATA 
SCAN  10,  60 


END  OF  ECL 
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b.  Alternative  File  Directive 


altfil 


ALTFIL 

This  directive  causes  the  scanner  to  read  from  the  "other" 

card  input  file.  For  convenience,  the  user  input  files  available  are  on 

SCOPE  files  TAPE4  and  TAPE5.  When  ALTFIL  is  encountered  in  the  ECL, 

I 4 I 

subsequent  input  is  read  from  TAPE  j 5 j depending  upon  which  file  is 
currently  being  read.  EXEMPT  starts  reading  TAPE5,  the  normal  system 
card  input.  When  the  first  ALTFIL  directive  is  encountered,  subsequent 
ECL  is  read  from  TAPE4.  When  an  END  or  ALTFIL  is  encountered  on  TAPE4, 
further  input  is  read  from  TAPES  until  another  ALTFIL  or  END  directive 
is  encountered.  All  input  is  transferred  to  TAPE3  before  processing. 
Note,  TAPE4  must  be  terminated  by  an  ALTFIL  directive. 

Example:  Assume  that  TAPE4  and  TAPES  are  structured  as  illus- 
trated in  (a)  and  (b)  of  the  illustration  shown  below. 


TAPES 

TAPE4 

TAPE3 

# 

• 

• 

• 

• 

ECL1 

ECL3 

ECL1 

# 

• 

• 

* 

• 

% 

ALTFIL 

ALTFIL 

ECL  3 

# 

# 

ECL2 

+ 

ECL4  

-►  ECL2 

• 

# 

• 

• 

• 

• 

ALTFIL 

ALTFIL 

ECL4 

• 

• 

• 

• 

• 

END 

DATA  2 

END 

• 

• 

. 

• 

• 

• 

DATA1 

END 

DATA1 

ALTFIL 

ALTFIL 

DATA2 

RUN 


c.  Run  Directive 

RUN 

This  directive  causes  the  LIP  to  terminate  processing  and 
execute  those  commands  encountered  from  the  beginning  of  the  ECL  or 
since  the  last  RUN  directive.  If  data  cards  are  to  be  read  using  the 
EXEMPT  READC  command,  they  must  be  placed  following  the  RUN  directive. 

Card  input  to  be  read  by  the  SLAM  or  user  provided  subroutines  appear  after 
the  7/8/9  card  following  the  ECL.  When  all  ECL  down  to  the  RUN  command 
has  been  executed,  control  returns  to  LIP  to  process  subsequent  commands. 
Note  that  a RUN  directive  followed  by  an  ALTFIL  directive  will  place  the 
contents  of  the  alternate  file  immediately  following  the  RUN  directive 
on  TAPE3,  the  input  file  used  for  execution.  This  is  illustrated  in  the 
following  example. 

TAPE5  TAPE4  TAPE3 

• • • 

• • • 

ECL1  DATA  ECL1 

• • • 

• • • 

RUN  4-  END  RUN 

ALTFIL  ALTFIL 

• • 

ECL2  DATA 

* • 

* • 

END  END 

ECL2 

END 

(a)  (b)  (c) 


End  Directive 


This  directive  is  used  to  separate  blocks  of  ECL  and  coded  data. 
When  encountered  by  the  LIP,  control  is  transferred  to  the  CLP  and  will  not 
be  returned  to  the  LIP.  As  illustrated  below,  DATA  would  be  read 
during  execution  of  ECL1.  Control  would  return  to  the  LIP  to  process  ECL2. 


r 


loop/label 


2.  COMMAND  LANGUAGE  PROCESSOR 

When  control  is  relinquished  by  the  LIP,  the  user's  ECL  has  been 
translated  into  a coded  execution  list.  The  data  reside  intact  on  TAPE3 
for  reading  during  execution  of  a READC  command  if  present.  There  are 
four  basic  types  of  operations.  These  are  control,  logical,  SLAM,  and 
arithmetic  operations. 

These  commands  will  be  treated  in  separate  sections  which  follow. 

a.  Command  LOOP /LABEL 


LOOP 


LABEL 


Name 

n 

Name 


n 


These  are  the  control  commands  present  in  this  release  of 
EXEMPT.  They  must  occur  in  pairs  and  cause  all  commands  between  them  to 
be  executed  the  number  of  times  specified  by  n.  The  name  or  number  used 
to  identify  the  LABEL  must  be  unique  (I.e.,  multiple  LOOPs  cannot  terminate 
on  the  same  LABEL).  However,  LOOPs  may  be  nested  to  any  level  as  long 
as  the  maximum  number  of  LOOP/LABEL  commands  does  not  exceed  100.  This 
is  the  internal  dimension  of  the  loop  control  table  and  may  be  easily 
changed  if  required.  Examples  of  both  simple  and  nested  loops  are 
presented  in  the  example  below.  Example  (a)  presents  a simple  loop 
using  an  integer  label  identifier  while  example  (b)  illustrates  nested 
loops  using  alpha  label  identifiers. 


ECL 

ECL 

LOOP  5,  10 

LOOP  OUTER, 

• 

LABEL  5 

• 

LOOP  NEST1, 

END 

LOOP  NEST 2, 

LABEL  NEST2 
LABEL  NEST1 
LABEL  OUTER 

(a) 

(b) 
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b.  EXEMPT  Logical  Connnands 

There  are  two  main  functions  fulfilled  by  the  logical  operations 
in  EXEMPT.  There  are  file  manipulation  and  data  set  I/O.  The  underlying 
commonality  of  these  operations  is  that  they  do  not  alter  the  internal 
data,  they  simply  format  them  and  read  them  into  EXEMPT  or  write  them  to 
some  device.  These  commands  will  be  further  classified  as  Random 
Access  File  Commands,  Sequential  File  Commands,  and  Data  Set  I/O  Commands. 

(1)  Random  Access  File  Commands 

EXEMPT  uses  FORTRAN  logical  unit  99  for  random  access 
I/O.  All  data  manipulated  internally  is  sorted  on  this  logical  unit. 

The  detailed  structure  of  this  logical  unit  is  presented  in  Section  III. 
The  user  may  have  up  to  20  files  identified  by  a name  or  number  on  this 
logical  unit.  Each  of  these  20  files  may  contain  up  to  54  data  arrays. 

A data  array  is  an  element  of  a data  set.  For  instance,  frequency  and 
time  domain  data  sets  contain  two  data  arrays,  the  dependent  (response) 
data  and  the  independent  (frequency  or  time)  data.  A typeless  data  set 
contains  only  one  data  set  (itself).  In  order  to  conserve  space,  whenever 
an  independent  data  array  is  stored  it  is  compared  with  existing  data 
arrays  of  the  same  type,  and  if  found  to  be  identical  to  an  existing 
array  it  will  not  be  stored.  This  prevents  clogging  up  the  file  with 
duplicate  copies  of  the  independent  data  which  are  usually  the  same  for 
many  data  sets. 

The  following  commands  are  to  aid  the  user  in  creating, 
maintaining,  and  modifying  his  random  access  data  base.  The  data  to 
reside  on  the  data  base  may  have  several  sources.  It  may  be  read  from 
cards,  generated  by  EXEMPT  itself,  or  interfaced  from  another  computer 
program  using  the  sequential  I/O  commands.  Usually,  some  combination  of 
these  methods  is  used.  When  EXEMPT  starts  execution,  it  checks  to  see 
if  a data  base  exists.  If  one  does,  it  is  loaded  into  EXEMPT.  The  user 
must  have  attached  the  data  base  as  SCOPE  local  file  TAPE99.  If  no  data 
base  exists,  EXEMPT  creates  one  for  use  during  this  particular  run.  If 
the  user  has  requested  that  TAPE  99  reside  on  a permanent  file  device 
(REQUEST,  TAPE99,  * PF.),  the  created  data  base  may  be  cataloged  as  a 
permanent  file  when  the  run  concludes  (CATALOG,  TAPE99 , MYFILENAME, 


FILE 


FILE 


(a)  File  Identifier 

I NAME  I 

I » I 

This  command  designates  the  random  access  file  to  be 
used  for  data  retrieval  and  storage.  All  data  generated  subsequent  to 
this  command  will  be  9tored  on  the  file  identified  by  NAME  or  N.  If 
data  already  exist  on  this  file,  they  will  be  made  available  to  EXEMPT. 

As  stated  earlier,  20  such  files  may  be  designated  by  the  user.  If  a 
file  is  never  identified,  the  file  name  TEMP  is  assumed. 

In  the  following  example,  assume  that  file  EXPDAT 
contains  experimental  data  to  be  used  and  it  is  desired  to  save  the  results 
obtained  on  a _‘ile  identified  as  F15ANL. 


FILE 


EXPDAT 


ECL  TO  RETRIEVE  DATA 


FILE  F15ANL 


ECL  TO  PERFORM  ANALYSIS 


$SUBSEQUENT  DATA  STORED 
$0N  FILE  NAMED  F15ANL 


END 
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SAVE 


(b)  Save  Data  Command 

SAVE  Name 


[• 


FILE 


\ Name  | "1 

| n | J 

This  command  will  result  in  the  data  associated  with 
the  name  specified  being  written  on  the  random  access  file  specified. 

If  no  file  is  specified,  the  current  file  is  used.  If  the  file  has  not 
been  previously  created  via  a FILE  command,  one  will  be  automatically 
creaced  before  the  data  are  transferred. 

This  command  is  generally  used  to  transfer  data  from 
the  working  file  TEMP  to  another  file  to  be  saved  on  a permanent  basis. 
All  data  sets  which  are  created  during  an  execution  will  be  stored  on. 
the  current  file.  Unless  specified  otherwise,  this  file  has  the  name 
TEMP  and  may  become  full  if  data  are  not  selectively  moved  to  another 
file  when  not  required  in  the  subsequent  commands.  This  is  illustrated 
below  where  the  file  TEMP  is  used  by  default  to  generate  data  set  ANT47A 
and  the  data  are  subsequently  transferred  to  file  ANTRSP  with  the  SAVE 
c ommand . 

ECL  TO  GENERATE  DATA  SET  AFTFRQ 


ANT47A  - INVFT  (ANTFRQ) 
SAVE  ANT47A,  FILE  - ANTRSP 


END 


The  data  set  name  (and  data)  ANT47A  will  no  longer 


PURGE 


(c)  Purge  Random  File  Command 


PURGE 


FILE 


| Name  ) 

I N I 


[ , SAVEFN] 


This  command  causes  all  data  stored  on  the  specified 
file  to  be  purged.  Unless  SAVEFN  is  specified,  the  file  name  will  also 
be  purged  from  the  master  file  directory.  In  practice,  this  command  can 
be  used  to  clean  up  a data  base  by  purging  data  no  longer  required.  It 
may  also  be  used  to  eliminate  the  data  of  file  TEMP  before  the  conclusion 
of  a run  to  prevent  subsequent  users  of  the  data  base  from  using  temporary 
data  saved  during  the  current  run. 
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READF 


(e)  Read  File  Data 

(1)  READF  Name  [.FILE  -lName)  ] 

L | n ) J 

(2)  READF  FILE  - j Name 

( n 

The  two  forma  of  this  command  perform  two  distinct 
functions  and  must  not  be  confused.  In  order  to  aid  in  the  discrimina- 
tion of  these  functions,  the  commands  will  be  discussed  separately. 

(1)  READF  Name  I", FILE  - j Name 

L I n 

This  command  causes  the  data  set  specified  to  be 
recovered  from  the  random  access  file  specified.  If  no  file  is  specified, 
then  the  current  file  is  searched  for  the  data  set.  A very  important 
action  takes  place  when  the  data  set  is  recovered.  First,  the  data  set 
is  stored  on  the  current  random  access  file  if  it  was  retrieved  from 
another  random  access  file.  Secondly,  the  global  independent  array 
corresponding  to  the  data  set  type  is  overwritten  with  the  independent 
data  of  the  data  set.  Thus,  if  the  data  set  is  a frequency  domain  type, 
FREQ  will  be  overwritten  with  the  frequency  values  associated  with  the 
data  set.  If  it  was  a time  domain  type,  TIME  would  be  overwritten.  The 
primary  purpose  of  this  command  is  to  allow  the  user  to  use  the  same 
FREQ  and  TIME  data  for  data  sets  which  reside  on  different  files.  The 
names  FREQ  and  TIME  may  not  be  used  directly.  This  is  due  to  the  fact 
that  these  are  special  names  and  there  may  be  several  arrays  on  each 
file  with  these  names.  EXEMPT  keeps  track  of  the  different  data  by 
linking  these  independent  arrays  to  a particular  response  data  array  to 
form  a complete  data  set. 

For  example: 

READF  ANT47A  FILE  = ANTRSP 

would  retrieve  the  data  associated  with  ANT47A  and  replace  the  global 
TIME  array  with  that  associated  with  data  set  ANT47A  on  file  ANTRSP.  If 
ANTRSP  is  the  current  file  name,  then  READF  ANT47A  would  simply  overwrite 
TIME  with  the  independent  data  of  ANT47A. 


The  second  form  of  this  command  is: 

(2)  READF  FILE  - J N®m®  | 

This  command  informs  EXEMPT  of  the  data  set  names  present 
on  the  file  specified.  During  LIP  execution,  all  data  set  names  are 
entered  in  a symbol  table;  however,  the  file  location  is  unknown.  When 
this  command  is  executed,  the  names  of  the  data  sets  on  the  specified 
file  are  recovered  and  if  a match  is  found  in  the  symbol  table,  EXEMPT 
Is  informed  that  the  data  associated  with  the  symbol  reside  on  the 
designated  file.  This  allows  the  user  to  have  data  sets  with  the  same 
name  on  different  files.  If  a data  set  is  referenced  before  it  is  defined 
either  by  appearance  in  a subroutine  call  or  as  the  result  of  an  operation, 
a fatal  error  will  occur.  This  is  analogous  to  the  use  of  an  undefined 
symbol  in  FORTRAN.  If  the  user  wishes  to  use  data  from  a previously 
generated  data  base,  a READF  command  for  each  of  the  files  containing 
the  data  sets  to  be  used  is  appropriate.  This  feature  also  allows  the 
use  of  different  data  referenced  by  the  same  name  at  different  points  in 
the  execution.  This  occurs  since  the  ECL  is  executed  sequentially  and 
the  READF  command  will  only  effect  those  occurrences  of  the  data  set 
name  which  follow  It.  This  is  illustrated  in  the  example  below. 

READF  FILE  - E4LIB 

ECL  SDATA  SETS  ON  E4LIB  WILL  BE  USED 

READF  FILE  - B1LIB 

ECL  $DATA  SETS  ON  B1LIB  WILL  BE  USED 

END 

NOTE:  Unique  data  set  names  contained  on  E4LIB  are  still  available  to  the 
ECL  following  the  READF  FILE-B1LIB  command.  However,  those  names 
that  are  the  same  will  be  available  only  from  B1LIB. 
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SAVRAN 


(f)  Save  Random  Access  Files 

SAVRAN  ^FILE  - | N*me  | j 

This  command  permits  the  random  access  data  base  to  be 
saved  on  a sequential  file  which  would  normally  be  a magnetic  tape.  This 
is  required  since  the  data  base  may  grow  to  a size  which  demands  a lot  of 
computer  resources.  Also,  this  provides  a backup  capability  for  the  data 
base.  If  the  file  is  specified,  it  must  be  present  as  a SCOPE  local 
file  with  the  name  specified.  If  the  numeric  designator  is  used,  TEMPn 
must  be  the  SCOPE  file  name.  If  n is  a single  digit  number,  it  must  be 
preceded  by  a zero.  If  the  file  designator  is  omitted,  the  default  name  is 
NEWDAT  and  a SCOPE  file  with  this  local  name  must  be  available. 

The  sequential  file  format  and  content  is  presented  as 

follows. 

Record  Contents 

1 Active  file  units  table  and  number  of  active  file  units 

2 File  header  and  director,  first  active  file 

3 Frequency/time  response  or  independent  array,  data 
set  1 

4 Frequency/ time  response  or  independent  array,  data 
set  2 


N+2 

N+3 

N+4 


Frequency/ time  response  or  independent  array,  data 
N (N  < 54) 

File  header  and  directory,  second  active  file 

Frequency /time  response  or  independent  array,  data 
set  1 


M+-N+3  Frequency/ time  response  or  independent  array,  data 

set  M (M  < 20) 
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GETRAN 


(g)  Recover  Random  Access  Data  Base 

1 Name  I 


GETRAN 


[• 


FILE  - 


, n 1 J 

This  command  allows  recovery  of  Che  random  access 
data  base  saved  by  the  SAVRAN  command.  The  same  rules  apply  for  the 
SCOPE  local  file  names,  except  the  default  name  is  OLDDAT.  GETRAN  can 
only  recover  files  when  they  are  written  exactly  as  specified  in  the 
SAVRAN  command. 


■ '7T/ 


LIST  INDEX 


(h)  List  File  Contents 


LIST  INDEX 


This  command  allows  the  user  to  list  the  file  names 
and  data  set  directories  for  the  random  access  data  base.  If  the  par- 
ticular file  is  designated  in  the  FILE  field,  the  index  is  restricted  to 
that  file  specified.  Note  that  the  contents  of  the  data  sets  are  not 
listed.  Another  LIST  format  is  presented  under  OUTPUT  commands  for  this 
feature.  The  output  from  this  form  identifies  the  file  name,  data  set 
name,  data  set  type,  location  of  independent  array,  and  the  data  set 
title. 
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must  contain  120  or  less  characters. 

If  the  Name  and  FILE  are  omitted,  the  title  is 
associated  with  the  current  file.  If  Name  and  FILE  are  both  specified, 
then  the  title  is  associated  with  the  data  set  on  the  specified  file. 

When  only  one  of  the  fields  is  missing,  the  title  will  be  associated  with 
the  field  provided. 


RENAME 


(2)  Data  Display  Commands 

Due  to  the  large  amounts  of  data  which  may  be  encountered 
In  an  execution,  output  Is  kept  to  a minimum  except  In  the  debug  mode. 
Almost  all  output  must  be  directed  by  the  user.  The  data  display  commands 
permit  the  user  to  obtain  a list  of  the  data  on  the  printer,  punched 
cards,  and  various  plot  devices.  These  actions  are  In  response  to  the 
LIST,  PUNCH,  and  PLOT  commands  discussed  In  the  following  sections. 


FORMAT 


\ 


> 


(a)  Data  Formatting 

( MP  ) 

FORMAT  Name  [ , Name ' Rl  l 

( DB  j 

This  command  directs  the  format  of  the  output  for 
the  data  sets  specified.  It  applies  only  to  complex  data  since  real 
data  are  rigidly  formatted.  The  command  will  specify  that  whenever  a 
data  set  is  displayed,  it  will  be  in  real/imaginary  (RI)  or  dB  (DB) 
format.  If  not  specified,  data  is  displayed  in  magnitude/phase  (MP) 
format.  If  RADS  is  specified,  the  phase  will  be  displayed  in  radian 
measure.  If  not  specified,  the  phase  will  be  in  degrees.  In  all  cases, 
phase  is  bounded  by  +180°. 


DGREES  | 
RAD  j 


LIST 


'i 


I 


(b)  List  Data  Set  Contents 

(1)  LIST  Name  1, 

(2)  LIST  CURFIL 

(3)  LIST  FILE  - 

Form  1 of  this  command  is  used  to  output  the  data 

associated  with  the  data  sets  Name  1,  Name  2 ...  The  data  will  be 

listed  sequentially  using  the  format  specified.  DB  implies  dB  and  the 
data  will  be  normalized  to  the  maximum  value  in  the  array.  RI  speci- 
fies listing  the  real  and  imaginary  components  separately  while  HP  directs 
the  output  to  be  in  magnitude  and  phase  format.  The  default  format  is 

MP.  These  formats  apply  only  to  complex  domain  data.  Noncomplex  data, 

being  real  values  only,  are  printed  as  they  exist  in  the  computer.  All 
typed  data  lists  will  contain  both  the  independent  and  dependent  data. 

Form  2 of  the  command  will  cause  the  data  associated 
with  all  data  sets  on  the  current  random  access  file  to  be  listed.  The 
output  format  used  will  be  that  specified  on  a FORMAT  directive.  If  no 
format  directive  is  specified,  the  data  will  be  displayed  in  magnitude/ 
phase  for  complex  data  and  real  for  noncomplex  data. 

Form  3 lists  the  same  data  as  form  2 for  the  file 


specified. 


PUNCH 


(c)  Punch  Data  Output 

PUNCH  Name  1 [.Name  2,...]  [.SINGLE]  [ , INDPNT ] 

This  command  will  punch  the  data  associated  with  the 
data  set  names  specified  In  all  cases,  the  data  from  Name  1 will  be 
punched  before  any  data  from  Name  2.  If  SINGLE  is  specified,  a single 
data  entry  will  be  punched  per  card.  If  SINGLE  is  not  specified,  there 
will  be  six  fields  per  card  containing  data  determined  by  the  data  set 
type  and  the  presence  of  the  INDPNT  directive.  If  INDPNT  is  specified, 
the  value  of  the  independent  variable  will  be  punched  with  the  value  of 
the  dependent  data.  If  INDPNT  is  not  specified,  only  the  response  data 
will  be  punched.  All  data  are  punched  in  12  column  fields  with  the  first 
column  of  each  field  left  blank.  The  data  are  punched  on  a 1XG11.5  format 
for  each  field  and  columns  76  through  80  will  contain  the  Julian  date. 

If  INDPNT  is  specified,  the  independent  data  will  be 
punched  preceding  each  value  of  the  dependent  data.  For  a frequency 
domain  data  set,  this  corresponds  to  three  numbers  for  each  element  of 
the  data  set  while  two  numbers  suffice  for  time  domain  data.  If  SINGLE  is 
also  specified,  this  is  the  limit  for  each  card.  However,  if  SINGLE  is 
not  specified,  two  elements  of  a frequency  domain  data  set,  or  three 
elements  of  a time  domain  data  set  will  be  punched  on  each  card. 

If  INDPNT  is  not  specified,  only  the  dependent  data 
will  be  punched.  This  corresponds  to  two  numbers  for  frequency  domain 
data  and  one  number  for  time  domain  or  typeless  data.  In  this  case,  if 
SINGLE  is  not  specified,  there  will  be  three  elements  per  card  for  frequency 
domain  data  and  six  elements  per  card  for  time  domain  or  typeless  data. 

Punched  header  cards  will  precede  each  punched  data 
set.  The  header  cards  will  contain  the  data  set  name,  number  of  points, 
and  the  data  set  title. 
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PLOT 

PRTPLT 


(d)  Plot  Display  Commands 


, Name , . . . ] 


* Hollerith  String  *]  [ ,TITLEY 


ILINLIN  \ 
L INLOG  I 
LOGLIN  I 
LOGLOG  ( 
DBLIN  | 
DBLOG  J 
* Hollerith  String  *] 


. I L.UUL1N  i 

) | Name  1 [ , Name  2 , . . . , Name , . . . ] [ , j LOGLOG  J ] 

|!>RIPLIi  DBLIN 

^ DBLOG  J 

(.TITLEX  = * Hollerith  String  *]  [,TITLEY  = * Hollerith  String  *] 

[ , OVRLAY ] 

These  commands  will  plot  the  data  sets  specified  on  the 
device  indicated.  PLOT  indicates  a METALIB  plot  file  is  to  be  created; 
whereas  PRTPLT  indicates  an  EXEMPT  printer  plot  is  to  be  placed  on  the 
output  file.  The  disposition  of  the  METALIB  plot  file  is  controlled  by 
the  user  prior  to  EXEMPT  execution.  The  plot  file  will  reside  on  LFN 
PI.TFIL  which  may  be  disposed  of  immediately  or  cataloged  as  a permanent 
file  for  later  disposition.  (See  Section  V for  an  example  of  the  METALIB 
file  disposition). 

To  receive  any  plots  through  the  METALIB  system,  the  user 
should  do  the  following  steps  prior  to  execution: 

(1)  Request  PLTPIL  to  be  a Q file. 

(2)  Dispose  PLTFIL  to  the  proper  plot  device. 

For  example,  to  receive  a CALCOMP  plot,  the  following  JCL 
cards  are  inserted  prior  to  EXEMPT  execution: 


REQUEST (PLTFIL, *Q) 
DISPOSE (PLTFIL , *MF»PCF) 
LGO. 


DISPOSE (PLTFIL , *MF=PCF , ST=ANY ) 
LGO. 


The  following  plot  devices  are  available  to  the  user: 


CODE 

DEVICE 

CHARACTERISTICS 

PCF 

CALCOMP  DRUM 

FULL 

SIZE,  .01  INCREMENTS, 

PLAIN  WHITE  PAPER 

PCG 

CALCOMP  DRUM 

FULL 

SIZE,  .01  INCREMENTS, 

RED  RULES  PAPER 

PG4 

GOULD  4800 

PVT 

VERSATEC 

PB1 

FR80 

B/W 

16  MM 

FILM 

PB3 

FR80 

B/W 

35  MM 

FILM 

PF8 

FR80 

48  X 

MICROFICHE 

i 


Also,  If  the  user  wishes  to  have  a printer  plot  of 
the  META  PLOT  file,  the  following  card  is  needed  after  EXEMPT  execution: 

DIRECT (PLTFIL, .PRINTER) 

When  PLTFIL  is  directed  to  the  printer,  the  META 
system  converts  the  file  for  printing  to  8 lines  per  inch  rather  than  the 
normal  6 lines  per  inch. 

The  axis  scaling  will  be  determined  by  the  LINLIN 
through  DBLOG  parameters.  The  first  acronym  refers  to  the  Y axis  and  the 
second  refers  to  the  X axis.  When  DB  is  specified,  the  data  will  be 
normalized  to  the  maximum  value  of  all  data  sets  in  the  list.  The  value 

plotted  will  be  20  log^  — ^ — . 

max 

The  axis  titles  are  controlled  by  the  TITLEX  and 
TITLEY  fields.  If  these  are  absent,  the  X axis  title  will  be  FREQUENCY 
or  TIME,  as  appropriate.  The  Y axis  title  will  be  DB  MAGNITUDE  on  LN 
MAGNITUDE,  or  MAGNITUDE  as  directed  by  the  axis  scaling.  If  the  TITLEX 
or  TITLEY  fields  are  provided,  the  Hollerith  string  between  the  asterisks 
will  appear  as  axis  titles.  Maximum  title  length  for  PRTPLT  is  30 
characters  and  for  PLOT  it  is  20  characters. 

Overplotting  of  data  sets  is  controlled  by  the  OVERLAY 
field.  For  PRTPLT,  only  one  data  set  per  plot  will  result;  however,  all 
plots  will  have  the  same  scale  for  easy  comparison.  For  PLOT,  up  to  four 
data  sets  may  be  overplotted  per  plot  frame.  The  reason  for  the  limitation 
is  that  for  CALCOMP  plots  there  is  only  enough  room  on  a frame  for  a plot  plus 
four  data  sets'  characteristics.  These  characteristics  include  the  data  set 
name  and  a 120  character  title.  More  than  four  may  be  overplotted;  however, 
only  the  first  four  data  sets  will  have  any  of  its  characteristics  printed 
on  the  plot  frame. 
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(3)  Sequential  File  Commands 

In  addition  to  the  random  access  file  capability,  EXEMPT 
also  will  manipulate  sequential  (typically  magnetic  tape)  files.  This 
type  of  file  is  the  primary  interface  with  other  user  programs.  The 
EXEMPT  subroutine  EXINTO  may  be  incorporated  in  a user  program  to  write 
an  EXEMPT  compatible  sequential  file.  This  routine  is  discussed  further 
in  Section  III.  In  general,  the  initial  input  data  for  EXEMPT  will  be 
experimental  or  generated  by  another  modeling  program.  The  user  may 
smooth  or  modify  the  data  before  they  are  written  to  tape  for  future 
incorporation  into  the  EXEMPT  data  base.  There  are  two  levels  of  control 
allowable  on  the  sequential  file.  Entire  data  sets  (header,  dependent, 
and  independent  data)  are  referenced  as  blocks.  Groups  of  blocks  are 
delineated  by  FORTRAN  end  of  file  marks  and  the  file  may  be  manipulated 
in  either  blocks  or  files.  This  permits  the  user  to  minimize  his  magnetic 
tape  library  since  he  can  direct  EXEMPT  to  position  the  tape  at  the 
correct  data  set  (block)  and  subsequently  read  the  data  and  associate 
them  with  a specific  internal  EXEMPT  symbolic  name. 

The  following  pages  present  a description  of  the  sequential 
I/O  commands  and  illustrations  of  their  use.  Detailed  illustrations  of 
the  sequential  file  structure  are  presented  in  Section  III. 
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FRQFIL 

TIMFIL 


(a)  File  Designation  Command 


FRQFIL  = 

( Name 
n 

TIMFIL  - 

j Name 
I n 

This  command  establishes  the  sequential  files  for 
frequency  and  time  domain  data.  These  data  must  be  stored  on  separate 
sequential  files.  All  subsequent  sequential  file  identifiers  may  then 
be  FRQFIL  or  TIMFIL  and  the  proper  user  file  will  be  referenced.  There- 
fore, in  the  following  command  descriptions,  FRQFIL  or  TIMFIL  may  be 
used  as  the  file  identifier.  The  file  identified  in  this  command  must 
be  attached  as  a SCOPE  local  file  with  the  same  name  identifier  as  the 
TEMPn  designator  when  n is  used.  When  n is  a single  digit,  it  must  be 
preceded  by  zero. 
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SKIPF 


(b)  Skip  File  Command 

SKIPF  jN*me|  [.  ±n] 

This  command  will  skip  the  number  of  FORTRAN  files 

determined  by  |n|  . If  N is  positive,  the  file  will  be  positioned  forward 

while  if  N is  negative,  the  file  will  be  positioned  backwards  to  a 

preceding  file.  If  N is  omitted,  it  will  default  to  +1.  The  file 

identifier  j darnel  must  be  attached  by  the  user  as  a SCOPE  local  file 
I n ) 

name.  If  the  numeric  file  identifier  n is  used,  the  local  file  name  is 
TEMPn.  If  Name  is  used,  then  the  local  file  must  have  the  same  name. 

For  example,  suppose  that  a user  has  two  sequential 
files  on  disk  for  convenience.  Then  the  SCOPE  control  cards  and  ECL 
would  resemble 

ATTACH,  NAMFIL,  SEQFIL1,  ID  - USRID. 

ATTACH,  TEMP23,  SEQFIL2,  ID  - USRID. 

EXEMPT. 


ECL 

SKIPF  NAMFIL  5 $ SKIP  F0RWARD  5 FILES 

SKIPF  23  1 $ SKIP  FORWARD  1 FILE 

ECL 

SKIPF  NAMFIL  -2  $ SKIP  BACK  2 FILES 
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SKIPB 


(c)  Skip  Blocks  Command 


SKIPB  j Name  j 

I n I 


| frqfil| 

I TIMFIL  j 


This  command  will  cause  EXEMPT  to  reposition  the  file 
by  N data  sets.  Positive  N results  in  forward  movement,  while  negative 
N will  reposition  the  file  to  a previous  data  set.  The  file  identifier 
may  be  a name  or  number  n and  must  be  a SCOPE  local  file  name.  If  the 
file  identifier  is  not  FRQFIL  or  TIMFIL,  then  one  of  these  must  be 
specified  to  Inform  EXEMPT  of  the  type  of  data  contained  on  the  file. 

Only  one  type  of  data  is  allowed  when  using  the  SKIPB  command. 

For  example,  the  following  SCOPE  and  ECL  would 
position  the  file  before  the  third  data  set  of  the  second  file  on  SCOPE 
local  file  FRQRSP . 


ATTACH,  FRQRSP,  USRPFNAME,  ID  - USERID. 
EXEMPT. 


FRQFIL  - FRQRSP 
SKIPF,  FRQFIL,  2 
SKIPB,  FRQFIL,  2 


If  the  FRQFIL  command  had  not  been  used,  then  the  ECL  would  be 


SKIPF,  FRQRSP,  2 
SKIPB,  FRQRSP,  2,  FRQFIL 
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READ 


(d)  Read  Sequential  File 

READ  Name  1 [,  Name  2,  ...  Name],  FILE  - J N*“e  J 

This  command  will  cause  the  number  of  data  sets  named 

in  the  list  to  be  read  from  the  designated  file  and  stored  on  the  current 
random  access  file  with  the  specified  names.  The  data  sets  will  be  read 

sequentially;  that  is,  the  first  will  be  identified  as  Name  1,  the  second 

as  Name  2,  etc.  The  type  of  data  and  format  are  specified  in  the  header 
record.  Subsequent  reference  to  the  data  sets  will  use  the  data  read 
from  the  sequential  file. 

Continuing  the  sample  from  the  SKIPB  command,  we 
could  read  the  next  three  data  sets  into  EXEMPT  with  the  ECL  instruction 
READ  RSP1,  RSP2,  RSP3  FILE  - FRQFIL 
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(e)  Read  Cards  Command 


READC 


R£ADC  Name  (,  Nl],  Name  [,  N2]  ... 

READC  Name , Name , ...  [ , N ] 

This  command  will  result  in  a coded  card  input 
read.  The  data  will  be  decoded  for  entry  into  the  specified  array 
storage  area.  If  the  number  N is  specified  with  the  array,  the  arrays 
are  filled  sequentially.  If  the  number  is  specified  at  the  end,  the 

arrays  are  filled  in  parallel.  The  equivalent  FORTRAN  statement  is 
shown  below. 


EXEMPT: 

READC  A,  5,  C,  8 

FORTRAN: 

READ*,  (A(I),  I - 1,  5), 

(C(I),  I * I,  8) 

EXEMPT: 

READC  A,  B,  C,  10 

FORTRAN: 

READ*.  (A(I) , B(I) , C(I) , 

1=1,  10) 

If  the  named  data  sets  have  independents,  then  those 
independents  are  read  as  well  as  the  imaginary  part  of  a frequency  domain 
data  set. 

The  cards  to  be  read  are  imbedded  in  the  ECL  and  on 
the  card  input  file  on  the  ALTFIL  file.  This  is  illustrated  below. 


READC  A,  B,  C,  10 


END 


DATA  FOR  A,  B,  C in  free  field  format 


END 
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WRITE 


(f)  Write  Sequential  File  ( J 

WRITE  Name  1 (.  Name  2,  ...  Name],  FILE  ” | N j 

This  command  will  write  the  data  associated  with  the 

data  sets  named  on  the  file  specified  in  EXEMPT  format.  This  format  is 

discussed  with  the  EXINTO  subroutine  in  Section  III. 

Continuing  the  example  from  the  READ  command,  assume 
that  RSP1,  RSP2,  and  RSP3  had  been  modified  by  the  execution  and  have  the 
same  names.  It  is  desired  to  save  the  new  data  in  place  of  the  old  data, 

and  they  can  be  the  last  data  on  the  file.  The  ECL  would  be 

SKIPB  FRQFIL  -3 
WRITE  RSP1,  RSP2,  RSP3 
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(g)  FORTRAN  End  of  File 


This  connand  will  write  a FORTRAN  end  of  file  (EOF)  on 
the  me  specified,  This  he  „3ed  t0  separate  data  ^ ^ 

° * Pesfflonlng  the  tape  ,1th  SCOPE  control  cards.  It  also  aids 

EXEMPT  in  positioning  the  file  efficiently  .hen  „r.  than  one  file  1, 
on  the  file. 


( 


• i 

Md 


rewind 


(h)  Rewind  Sequential  File 

REWIND  J N*me  | 

This  command  will  rewind  the  designated  file  and 
reposition  it  at  the  beginning  of  the  first  file.  The  FORTRAN  rules 
apply  in  that  if  the  last  operation  on  the  file  was  a WRITE,  an  EOF  will 
be  written  before  the  file  is  rewound.  If  the  last  operation  was  a 
READ,  the  file  is  simply  repositioned. 


(i)  Return  Sequential  File 
DROP  JN“ej 

This  command  returns  the  designated  file  to  the 
system.  Its  primary  purpose  is  to  prevent  EXEMPT  from  tying  up  a tape 
servo  after  the  data  have  been  read.  The  DROP  command  will  free  the 
unit  for  other  users  while  EXEMPT  executes. 


I 
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(4)  EXEMPT  Library  Functions 

EXEMPT  contains  functions  to  perforin  certain  numeric 
operations  on  previously  defined  data  and  to  generate  analytical  data 
sets  to  support  linear  systems  analysis.  All  of  these  functions  may 
apppear  in  arithmetic  expression.  When  used  with  an  arithmetic  expression, 
the  total  function  text  must  be  enclosed  within  parentheses.  The  two 
forms  are  illustrated  below. 

Name  » FUNCTION  (Name  1,  Name  2),  PI  ■ n,  P2  ■ m 
Name  - A * (FUNCTION  (Name  1,  Name  2),  PI  - n,  P2  - m) 

In  the  first  example,  the  function  will  generate  the  data 
set  name  by  operating  data  sets  Name  1 and  Name  2 using  the  parameters  PI 
and  P2.  In  the  second  example,  the  data  generated  by  the  function  will  be 
multiplied  by  A before  being  stored  as  Name. 

As  a rule,  the  parameters  in  the  list  are  optional  and  if 
no  parameters  are  associated  with  the  command  format,  the  command  need 
not  have  the  enclosing  parenthesis  when  used  in  an  arithmetic  expression. 

If  parameters  are  in  the  command  format  but  not  used,  the  enclosing  parenthesis 
are  required. 


k 
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(a)  Data  Interpolation 

Name  - INTERP  (Name  1 [,  Name  2])  [,  NPTS  ^ n]  [,  PT1  - v,]  [PT2  - v] 

This  command  will  Interpolate  the  dependent 
data  associated  with  Name  1 onto  a new  independent  data  base  associated 
with  Name  2 and  store  the  result  as  Name.  If  Name  2 is  not  provided, 
then  the  global  frequency  or  time  array  will  be  used  depending  on  the 
type  of  data  in  Name  1.  The  interpolation  is  an  Aitken-Lagrange  polynomial 
of  degree  specified  NPTS  * n.  NPTS  may  have  any  value  from  1 to  10.  The 
default  value  is  NPTS  ■ 1,  or  linear  interpolation.  The  parameters  PTl 
and  PT2  specify  the  first  and  last  values  of  the  new  independent  data. 

If  PTl  is  not  provided,  interpolation  will  start  with  the  first  point  of 
the  new  independent  data.  If  PT2  is  not  provided,  interpolation  will 
continue  to  the  last  point  of  the  new  independent  data. 

If  the  new  independent  data  of  Name  2 extend 
beyond  the  domain  of  the  data  in  Name  1,  a warning  will  be  given  to  the  user. 
The  new  value  will  be  the  last  value  in  the  domain  of  Name  1.  That  is,  the 
extended  data  set  will  contain  constant  values  in  the  extended  region. 

In  the  following  example,  assume  VOLTS  is  a frequency 

domain  response  function  for  unequally  spaced  values  of  frequency  from  1 
to  100  MHz.  If  a Fast  Fourier  Transform  (FFT)  is  to  be  used,  then  the 

independent  data  must  be  equally  spaced  and  start  at  a value  of  zero. 


ECL 

FREQ  = 0. , LIN512,  100E6 
VOLTS  = INTERP  (VOLTS) 

END. 

The  resulting  values  associated  with  VOLTS  for 
frequencies  between  0 and  1 MHz  will  be  the  original  value  at  1 MHz. 
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(b)  Data  Editing 


EDIT 


1 

1 


i 

i 

; 

; 

! 


f 


Name  - EDIT  (Name  1)  [,  PT1  ■ v]  [,  PT2  - v] 

This  command  removes  bad  and  redundant  data  from  the 

data  set  Name  1 and  stores  them  as  data  set  Name  2.  The  parameters  PT1  and 

PT2  have  the  same  meaning  as  the  INTERP  function. 

ZERO  is  a user  definable  parameter  intended  to  represent 

a reasonable  criteria  for  evaluating  equality  between  floating  point  numbers. 

ZERO  default  value  is  l.E-10. 

Redundant  points  occur  when  two  consecutive  values 

of  the  independent  data  I and  I are  such  that 

n n+l 


< ZERO 


Bad  data  points  occur  when  two  consecutive  values 
of  the  dependent  data  D^  and  Dn+^  are  such  that 


> — — 

ZERO 
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MERGE 


(c)  Merging  Data  Sets 
Name  “ MERGE  (Name  1,  Name  2) 


i MAXVAL  | 
j MINVAL  j 


This  function  allows  the  user  to  combine  data  from 
two  data  sets  into  one  data  set.  The  source  data  associated  with  Name  1 
and  Name  2 will  be  combined  in  ascending  order  of  their  independent  variable. 
When  two  values  of  the  independent  variable  are  found  to  be  the  same  (within 
an  error  of  ZERO),  the  action  directed  by  MAXVAL,  MINVAL,  or  the  absence 
of  these  parameters  will  be  taken.  If  MAXVAL  is  specified,  the  maximum 
dependent  value  will  be  chosen.  If  MINVAL  is  specified,  the  minimum  will  be 
chosen.  If  neither  parameter  is  present,  the  resultant  value  will  be  the 
average  of  the  values. 

From  the  INTERP  example,  VOLTS  could  be  extended  to 
zero  frequency  using  a MERGE  function.  Suppose  NULL  is  a data  set  containing 
a zero  response  for  frequencies  of  0 and  1 MHz.  Then  the  following  ECL 
would  generate  VOLTS  as  an  equally  spaced  data  set  between  0 and  100  MHz. 


FREQ  - 1. , LIN101,  100. E6 
VOLTS  - INTERP  (VOLTS) 

VOLTS  - MERGE  (VOLTS,  NULL) 
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FT 

INVFT 

(d)  Fourier  Transform 

Name  * FT  (Name  1 [,  Name  2]) 

Name  - INVFT  (Name  1 [,  Name  2]) 

The  commands  will  transform  the  data  sets  Name  1 to 

the  time  (INVFT)  or  frequency  (FT)  domain.  If  Name  2 is  provided,  the 

output  data  set.  Name,  will  have  the  same  independent  data  as  Name  2.  If 

Name  2 is  not  provided,  the  output  independent  data  will  be  generated  using 

the  Nyquist  criteria  on  the  independent  data  of  Name  1.  For  example,  if 

I and  I . are  the  extremes  of  the  independent  data  associated  with  Name  1 
max  min  r 

which  contains  N data  points,  then  the  independent  data  of  Name  will  have 
increments  6 computed  from 


I 

! 


Name  1 are: 


, 1 _ 

2 ft I - I . )/N  \ 

( max  min  I 

If  Name  2 is  not  provided  and  the  independent  data  for 
1 - Equally  spaced 


then  the  output  data  Name  will  be  generated  using  the  FFT.  This  is  the  only 
method  to  employ  the  FFT.  All  other  input  will  cause  the  use  of  a full  Fourier 
integration.  The  INTERP  function  is  provided  to  allow  the  user  to  transform 
data  from  unequally  spaced  independent  arrays  to  equally  spaced  independent 
arrays. 

In  this  release  of  EXEMPT,  there  is  no  checking 
performed  to  assess  the  validity  of  the  transform. 

In  the  following  ECL  examples,  EQUAL  is  assumed  to 
be  an  equally  spaced  frequency  domain  data  set  with  a zero  frequency  value. 

UNEQ  is  an  unequally  spaced  time  domain  data  set. 


EQTIM  = INVFT  (EQUAL) 

NTPUNQ  = INTERP  (UNEQ,  TIME) 
EQFRQ  « FT  (NTPUNQ) 

FRQUNQ  = FT  (UNEQ) 


I 


I 

! 
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EQTIM  will  be  an  equally  spaced  time  domain  data  set 
generated  by  the  FFT.  NTPUNQ  will  be  the  data  from  UNEQ  interpolated 
onto  the  global  TIME  array  which  we  will  assume  has  a zero  value  and  is 
equally  spaced.  Consequently,  the  FFT  will  generate  EQFRQ  as  an  equally 
spaced  frequency  domain  data  set  with  a zero  frequency  value.  FRQUNQ  will 
be  generated  using  the  full  Fourier  transform  of  the  data  in  UNEQ. 


TRDEXP 

(e)  Double  Exponential  Transform 
Name  - TRDEXP  (VI,  V2,  V3) 

This  function  generates  the  Fourier  transform  of: 


for  the  frequency  values  in  the  global  frequency  array  FREQ. 

The  output  data  set  Name  is  a frequency  domain  data 
set  which  may  be  used  as  any  other  frequency  domain  data  sets.  The  primary 
purpose  of  this  function  is  to  allow  the  user  to  have  an  analytic  function 
commonly  encountered  in  EMP  or  EMC  analysis.  With  the  proper  values  of  V2 
and  V3,  this  function  can  represent  either  an  EMP  or  lightning  pulse. 


(f)  Damped  Sine  Transform 


TRDSIN 


Name  - TRDSIN  (VI,  V2,  V3) 

This  function  generates  the  Fourier  transform  of: 
-V2t 

VI  e c sin  (2irV3t) 

for  the  frequency  values  in  the  global  frequency  array  FREQ. 

The  output  data  set  Name  is  a frequency  domain  data 
set.  A damped  sine  wave  is  a commonly  encountered  excitation  in  EMP  and 
EMC  analysis.  The  output  data  may  be  used  as  an  analytic  excitation  for 
system  analysis. 


(g)  General  Polynomial  Functions 


Name  - POLY  (VI,  V2,  V3,  V4) 

This  function  Is  used  to  generate  a complex  frequency 
dependent  function  which  can  represent  a factor  of  a polynomial  or  two 
terms  of  a polynomial.  POLY  generates  a real  and  imaginary  part  with  user 
designated  factor  and  power  of  to  (frequency  in  radians  per  second).  This 
method  will  calculate  factors  which  can  be  multiplied  to  form  a transfer 
function.  POLY  can  also  be  used  for  unfactored  polynomials  by  pairing 
real  and  imaginary  terms  of  to.  Each  pair  can  produce  a term  which  can  be 
added  to  form  a polynomial,  to  is  the  radian  frequency  derived  from  the 
global  frequency  array  as  2ir*FR£Q. 

V 2 V4 

POLY (VI,  V2,  V3,  V4)  - Vico  + j V3co 

Example: 

It  is  standard  notation  to  represent  a transfer  function  by: 


1 

(s  + a)  (s  + b) 

The  transfer  function  can  be  written  and  stored  in  TR  by: 

AB  - A*B 
C - A + B 

TR  - 1. / (AB-POLY (l.,2.,0.,0. )+POLY (0.,0.,1.,0.)*C) 
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(h)  Phase  Shift  Function 

Name  - SHIFT  (Namel,  t) 

This  function  is  used  to  shift  the  phase  of  the 
frequency  domain  data  set  Name  1 by  performing  the  following  operation 

. Name  » Namel,  * e^wt 


SHIFT 


for  those  values  of  w associated  with  the  frequency  data  of  Namel.  If 

f(t)  and  F(w)  are  transform  pairs,  then  f(t  + t ) and  F (w)  e^wto  are 

o 

also  transform  pairs.  This,  the  SHIFT  function  allows  the  user  to  adjust 
phase  without  recomputing  the  transform. 
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(i)  Data  Set  Typing 


TYPE  Name 


1 TIME 


| FREQ  \ 

This  declarative  statement  will  explicitly  define  a 
data  set  identified  by  Name  to  be  a frequency  or  time  domain  function. 

If  a data  set  is  not  referenced  in  a type  declaration,  it  will  be  typed 
implicitly.  Implicit  typing  occurs  when  a data  set  name  is  referenced 
as  the  resultant  of  a function  or  arithmetic  operation.  In  an  arithmetic 
operations,  the  first  independent  data  set  encountered  will  be  associated 
with  the  resultant.  If  none  are  encountered  and  the  resultant  is  not 
typed,  no  type  will  be  associated  with  it;  however,  if  none  are  encountered 
and  the  resultant  is  typed,  the  appropriate  global  independent  data  will 
be  associated  with  it. 

Attempting  to  perform  an  operation  or  evaluate  a 
function  with  improperly  typed  data  will  result  in  an  EXEMPT  fatal  error. 


. 


(j)  Data  Initialization 


Keyword 


Value  1 


, Value  2,  Value  3 


This  command  can  be  used  to  load  data  into  an  indepen- 
dent or  typeless  data  set.  The  most  common  use  will  be  to  initialize 
the  global  TIME  or  FREQ  data.  These  data  sets  may  be  redefined  anytime 
during  execution.  The  data  starts  at  value  1 and  increments  according 
to  the  specification.  The  meaning  of  the  specifications  are: 

LINn  Load  data  from  the  first  value  to  the  second 

value  with  n linear  steps 

LOGn  Same  as  LINn  but  with  logarithmic  steps 

REPn  Repeat  the  previous  value  n times 

NXTn  Repeat  the  next  value  n times 

When  initializing  the  arrays  TIME  and  FREQ,  the  first 
and  last  values  need  not  be  specified  for  the  LINn  and  LOGn  option.  As  can 
be  seen  in  Table  3,  minimum  and  maximum  values  for  TIME  (TIMMIN  and  TIMMAX) 
and  FREQ  (FRQMIN  and  FRQMAX)  are  already  defined.  Additionally,  these 
values  as  well  as  all  other  default  valued  keywords  may  be  redefined  by  a 
simple  replacement  expression  (i.e.,  A = value). 

Examples: 

FREQ  =*  LIN50 

This  command  loads  the  array  FREQ  with  50  equal 

linear  steps.  The  minimum  value  (FMIN)  is  0.  and  the  maximum  value 

(FMAX)  is  30.E+6. 

The  command: 

FREQ  - .01,  LIN10,  .1,  LIN10,  1.,  LIN10,  10.,  LIN10,  100. 
will  load  the  following  data  into  the  global  frequency  array. . 

.01,  .02,  .03,  .04,  .05,  .06,  .07,  .08,  .09,  .1, 

.2,  .3,  .4,  .5,  .6,  .7,  .8,  .9,  1., 

2.,  3.,  4.,  5.,  6.,  7.,  8.,  9.,  10., 

20.,  30.,  40.,  50.,  60.,  70.,  80.,  90.,  100. 
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Table  3 

KEYWORD  DEFAULT  LIST 


KEYWORD 

NAME 

INITIAL 

VALUE 

FUNCTION 

FREQMAX 

30.E+6 

MAXIMUM  FREQUENCY  DATA  VALUE 

FREQMIN 

0. 

MINIMUM  FREQUENCY  DATA  VALUE 

TIMMAX 

l.E-5 

MAXIMUM  TIME  DATA  VALUE 

TIMMIN 

0. 

MINIMUM  TIME  DATA  VALUE 

ZERO 

l.E-10 

CLOSENESS  TEST  FOR  TIME  OR 
FREQUENCY  DATA 

FRQFIL 

CONTAINS  THE  FREQUENCY  ONLY 
SEQUENTIAL  FILE  UNIT 

TIMFIL 

” 

CONTAINS  THE  TIME  ONLY 
SEQUENTIAL  FILE  UNIT 

PAGLIN 

51 

DEFAULT  VALUE  FOR  THE  PAGE 
LENGTH  IN  LINES  PER  PAGE 

PAGWID 

81 

DEFAULT  VALUE  FOR  THE  PAGE 
WIDTH  IN  COLUMNS  PER  PAGE 

R1-R20 

0 

GLOBAL  REAL  ARRAY  FOR  USER 

SUBROUTINES 
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(k)  Mathematical  Expressions 

Quantities  used  in  constructing  the  expression 
include  numerical  constants,  symbolic  constants,  function  references, 
data  sets,  and  global  variables.  Also  available  are  a variety  of  mathe- 
matical symbols  and  several  mathematical  functions. 

The  expression  is  written  using  constructions  similar 
to  those  used  in  FORTRAN.  EXEMPT  has  the  capability  to  evaluate  mathe- 
matical expressions  containing  data  set  names,  subroutine  names,  constants, 
and  arithmetic  functions. 

1)  Arithmetic  Operations 

The  arithmetic  symbols  +,  -,  *,  /,  and  **,  are 
used  to  denote  addition,  subtraction  or  negativity,  multiplication, 
division,  and  exponentiation,  respectively. 

All  arithmetic  operators  must  explicitly  appear 
in  mathematical  expressions.  No  operations  are  assumed  by  juxtaposition 
of  quantities. 

To  order  the  sequence  of  operations,  parentheses 
may  be  used.  For  example,  in  division,  only  the  quantity  immediately 
following  the  division  sign  is  considered  to  be  in  the  denominator. 

Compound  denominators  must  be  enclosed  in  parentheses.  These  examples 
illustrate  the  point. 


X6  = R1/.5*R3 
X7  - R1/(.5*R3) 


is  interpreted  as  X6 
is  interpreted  as  X7 


negative  of  a quantity: 


The  minus  sign  may  be  used  to  indicate  the 


-N(l) 


I 


j#  "7"' 


:l 


Exponentiation  is  indicated  by  the  symbol  **,  where 
the  quantity  appearing  immediately  after  the  symbol  is  the  exponent. 
Exponents  may  be  constant  or  variable.  Compound  exponents  must  be 
enclosed  in  parentheses.  Exponents  consisting  only  of  integer  constants 
should  be  written  without  the  decimal  point. 

Data  sets  contained  in  expressions  must  have  the 
same  associated  independent  data.  However,  a data  set  that  is  neither 
frequency  nor  time  dependent  may  be  mixed  within  an  expression. 

There  are  two  categories  of  mathematical  operations. 
If  a lower  case  letter  represents  a simple  variable  or  constant  (real 
or  complex),  and  upper  case  letter  represents  an  array  (real  or  complex), 
and  a ° represents  the  operators  +,  -,  *,  and  /,  then  the  categories  may 
be  described  as  follows. 

a.  First  type.  The  arguments  may  be  simple  variables,  arrays,  or 
both. 

c • a°b 
C - a°B 

■ (a#b^,  * •••»  a°bn) 

C - A°B 

- (a1°b1,  a2°b2 W 

b.  Second  type  (exponentiation).  The  arguments  may  be  simple 
variables  or  a simple  variable  and  an  array. 

c * a**b 
C - A**b 

■ (a  **b,  a **b,  ....  a **b) 

n 

2)  Parentheses 

Left  and  right  parentheses  are  available  for 
grouping  purposes.  Many  levels  of  parentheses  may  be  employed  to  a 
maximum  of  ten.  It  is  important  to  check  for  proper  pairing  of  left  and 
right  parentheses  so  that  ambiguous  or  incorrect  groupings  are  avoided. 

An  example  of  correct  usage  is: 

XI  - Rl/ (R2*(N+R3) ) 
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which  represents  the  expression 


XI 


R1 


R2 (N+R3) 

3)  Mathematical  Functions 

The  following  mathematical  functions  are  available 
for  use  in  mathematical  expressions: 


ABS(x) 

Absolute  value  of  x 

ATAN(x) 

Inverse  tangent  of  x,  result  ii 

ATAN2 (x,y) 

Inverse  tangent  of  y/x,  result 

CMPLX(x.y) 

Complex  function 

CONJG(x) 

Complex  conjugate  of  x 

COS(x) 

Cosine  of  x,  x in  radians 

EXP(x) 

Exponential  function  of  x 

ALOG(x) 

Natural  logarithm  of  x 

MAX(x.y) 

Maximum  value  of  x and  y 

MIN(x.y) 

Minimum  value  of  x and  y 

MOD(x,m) 

x modulo  m 

SIGN(x,y) 

Magnitude  of  x with  sign  of  y 

SIN(x) 

Sine  of  x,  x in  radians 

SQRT(x) 

Square  root  of  x,  x positive 

TAN(x) 

Tangent  of  x,  x in  radians 

Function  arguments  must  be  enclosed  in  parentheses. 
All  EXEMPT  functions  may  also  be  included  in  arithmetic  statements.  EXEMPT 
functions  with  parameter  list  must  be  totally  enclosed  in  parenthesis. 

Also,  if  a data  set  is  complex,  EXEMPT  will  perform  the  FORTRAN  complex 
function. 

There  are  three  categories  of  functional  opera- 
tions. If  a lower  case  letter  represents  a simple  variable  or  constant 
(real  or  complex)  and  an  upper  case  letter  represents  an  array  (real  or 
complex),  then  the  categories  may  be  described  as  follows: 

a.  First  type.  If  the  argument  is  a simple  variable,  the  result 

is  a simple  variable.  If  the  argument  is  an  array,  the  function 
is  applied  to  each  element  of  the  array  separately,  and  the 
result  is  an  array. 
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b. 


: | 

! 

. 


i 

\ 


C. 


b » f(a) 

B - f (A) 

- {f(a  ),  f(a  ),  ....  f (a  ) } 

i i-  n 

Second  type.  Functions  with  two  arguments.  If  both  arguments 
are  simple  variables,  the  result  is  a simple  variable.  If 
either  argument  is  an  array,  the  result  is  an  array, 
c - f (a,b) 

C - f (A,b) 

- { f (a^,b) , f(a2,b),  ....  f(an>b)> 
or 

C - f (A,B) 

■ {f(a1,b1),  f(a2,b2),  ....  f(an,bn)} 

Third  type.  The  argument  must  be  an  array  and  the  result  is 
an  array  or  both  arguments  must  be  arrays  and  the  result  is  an 
array. 

B = f (A) 

C - f (A,B) 


(5)  User  Subroutines 

EXEMPT  is  provided  with  the  capability  to  interface  with 
user  provided  software  using  two  methods.  The  first  method  involves 
replacing  existing  EXEMPT  subroutines  with  the  user  subroutine  of  the 
same  name.  There  are  45  of  these  available  and  they  are  identified  as 
subroutine  SUBA  SUBZ,  OPRO  ■+•  0PR9,  and  MYCOMO  -►  MYC0M9.  The  interface 
for  these  is  presented  in  Section  IV.  They  may  be  treated  as  any  other 
EXEMPT  function  with  respect  to  usage.  The  user  must  replace  the  EXEMPT 
subroutine  with  his  own  of  the  same  name  using  the  SCOPE  COPYL  command. 

The  identifiers  SUBa  and  OPRn  may  be  referenced  directly;  however,  the 
MYCOMn  is  a LIP  equivalency  statement.  The  user  defines  MYCOMn  as  a 
symbolic  name,  subsequent  reference  to  the  name  will  result  in  subroutine 
MYCOMn  being  called  by  EXEMPT.  This  is  illustrated  below.  When  the 
name  GTD  is  referenced  after  defining  MYC0M5  as  GTD,  MYC0M5  will  be  called 
to  perform  the  user's  operations  on  the  input  data. 

MYC0M5  = GTD 

ECL 

ANTRSP  = GTD  (FREQ) 

The  user  must  have  provided  a subroutine  MYC0M5  interfaced 
as  described  in  Section  IV. 

The  second  method  for  implementing  the  user's  software  is 
to  append  the  entry  point  for  a NAME  subroutine.  This  is  done  by  replac- 
ing the  EXEMPT  subroutine  CNGSUB  on  the  relocatable  binary  file  by  using 
the  COPYL  SCOPE  utility  with  the  append  option  specified  (see  example, 
Section  V).  The  replacement  must  include  an  EXTERNAL  declaration  followed 
by  a call  to  SETUP.  The  arguments  to  SETUP,  shown  below,  include  the 
user  subroutine  name  in  Hollerith  form  followed  by  the  subroutine  name. 

As  many  pairs  of  arguments  are  required  as  there  are  user-named  subroutines, 
up  to  a maximum  of  36  pairs.  \ 
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CALL 


SUBROUTINE  CNGSUB 

EXTERNAL  RUN1,  RUN 2 

CALL  SETUP (4HRUN1,RUN1,4HRUN2,RUN2) 

RETURN 

END 

. Further  information  and  detailed  description  of  the 
subroutine  uses  are  given  in  Section  IV.  It  is  recommended  that  the 
user  read  Section  IV  before  trying  a user  subroutine. 

The  general  types  of  calls  are  illustrated  below  where  A, 

B,  and  C represent  data  set  names  and  Name  is  used  as  a user  subroutine 
name. 


SECTION  III 

EXTERNAL  FILE  STRUCTURES 


There  are  two  data  file  structures  that  the  user  encounters.  These 
are  binary  sequential  files  and  random  access  files.  These  file  struc- 
tures are  used  automatically  through  the  command  execution.  Also,  it  is 
possible  for  a user  to  format  a binary  sequential  file  from  an  external 
program  for  U3e  on  EXEMPT. 

A master  index  will  contain  the  information  needed  for  accessing 
each  data  set  on  the  random  access  file.  The  master  index  points  to  a 
subindex  directory  for  the  particular  data  set. 

1.  EXEMPT  SEQUENTIAL  BINARY  FILE  FORMAT 

One  of  the  ways  in  which  the  EXEMPT  executive  computer  code  can 
exchange  data  with  other  independent  program  is  through  standard  format 
sequential  binary  files.  Since  the  files  are  binary  representations  of 
floating  point  numbers,  information  interchange  is  limited  to  machines  with 
the  same  word  length  and  floating  point  number  structure.  Currently, 

EXEMPT  standard  formats  are  defined  for  recording  frequency  domain  or 
time  domain  signals  along  with  some  title  information. 

The  external  unit  names  that  are  used  for  these  files  are  determined 
by  the  user.  The  unit  names  can  be  integer  numbers  between  1 and  99  or 
standard  FORTRAN  file  names  in  a 6L  format.  !f  the  files  are  referenced  by 
number,  the  external  name  for  these  files  will  be  TEMPXX.  The  files  for 
numbers  between  1 and  9 will  be  written  as  TEMP01  through  TEMP09. 

Therefore,  any  set  of  numbers  can  be  used  with  these  files  and  they  will 
not  be  confused  with  standard  FORTRAN  unit  numbers. 

The  structure  of  a sequential  file  is  in  terms  of  data  set.  Motion 
of  the  sequential  file  is  done  by  sets  rather  than  records.  A READ  will 
read  one  data  set  and  a WRITE  will  write  one  data  set. 

The  first  record  of  a set  will  contain  the  number  of  points,  data 
description,  date,  program  name,  and  the  title  information.  A time  domain 
data  type  will  have  two  other  records  and  a frequency  domain  data  type 
will  contain  three  other  records.  The  structures  of  these  data  types  are: 


Frequency  Domain  Data  (Four  Records) 

1.  Record  1:  17  words  long 

(a)  One  word  integer  (N)  indicating  number  of  words  in 
each  of  the  following  records. 

(b)  One  word  containing  data  description  in  CDC  display 


(c)  One  word  containing  data  description  of  data  values 
are  8HMAGPHASE,  8HREALIMAG,  7HDBPHASE. 

(d)  One  word  containing  a date  in  display  code. 

(e)  One  word  containing  display  code  or  name  of  program 
which  wrote  data  file  (e.g.,  7HSCEPTRE,  6HYSYNAP,  etc.). 

(f)  12  words  of  title  information  in  CDC  display  code. 

2.  Record  2:  N words  long  containing  the  array  of  frequency 
points. 

3.  Record  3:  N words  long  containing  the  real  part  or  ampli- 
tude of  a frequency  function. 

4.  Record  4:  N words  long  containing  the  imaginary  part  or 
phase  of  a frequency  function. 

b.  Time  Domain  Data  (Three  Records) 

1.  Record  1:  17  words  long;  same  as  for  frequency  data. 

2.  Record  2:  N words  long  containing  the  array  of  time  points. 

3.  Record  3:  N words  long  containing  the  values  of  the  time 
function. 

When  writing  an  external  file,  the  user  should  declare  his  file  to 
be  record  type,  RT=W  through  the  FILE  and  LDSET  SCOPE  utilities  or  use 
BUFFER  IN/OUT. 

To  properly  write  an  EXEMPT  formatted  sequential  file  as  described 
above,  a subroutine  EXINTO  is  available.  EXINTO  is  a FORTRAN  subroutine 
which  is  callable  by  other  FORTRAN  routines.  Its  purpose  is  to  write 
EXEMPT  format  data  files. 

The  proper  call  is: 


CALL  EXINTO  (ITYPE.LFN, YR, YI , FT, N, IOP, TITLE, IERR, NAME) 
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7>"—  : 


Arguments  are  as  follows: 

ITYPE:  A two-word  integer  array  identifying  the  type  of  data. 

ITYPE(l)  - 0 for  time  data  or  1 for  frequency  data 
ITYPE (2)  - 8HREALIMAG  or  8HMAGPHASE  or  7HDBPHASE 
as  appropriate  to  describe  frequency  data. 

LFN:  Logical  file  name  of  the  file  on  which  data  are  to  be 

written.  It  is  either  an  integer  from  1 to  99  or  a name 
with  a maximum  of  six  characters  in  L format.  The  system  file 
name  for  an  integer  specified  file  will  be  TEMPnn  where 
nn  is  the  integer.  The  file  name  must  be  declared  on 
the  FORTRAN  program  card.  To  reduce  memory  requirement, 
set  buffer  size  to  0. 

YR:  Data  array  containing  real  part  (or  magnitude)  of  a fre- 

quency dependent  function  or  amplitude  of  time  dependent 
function. 

YI:  Data  array  containing  imaginary  part  (or  phase)  of  a 

frequency  dependent  function  (not  used  by  time) . 

FT:  Data  array  containing  frequency  or  time  values. 

N:  Number  of  points  in  each  array. 

IOP:  Indicates  type  of  operation. 

1- Rewind 

2- Skip  data  block 

3- Backspace  data  block 

4- Skip  file 


5- Read  data  block 

6- Write  data  block 

7- Write  EOF 

8- Drop  File 

TITLE:  Up  to  12  words  of  Hollerith  information. 

IERR:  Return  error  codes: 

1- ITYPE  not  0 or  1 

2- Illegal  operation  code 

3- Illegal  file  name 

4- EOF  encountered 

5- End  of  information  (EOI)  encountered 
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6- Read  parity  error 

7- End  of  record  (EOR)  encountered,  data  array  not  N words  long 

8- Write  parity  error 

NAME:  Name  of  program  calling  EXINTO  in  left  justified  display 

code  (e.g.,  7HSCEPTRE  or  6HYSNAP) 

Example: 

CALL  EXINTO  (0, 3 , AMP, 0, TIME, 679, 6, TITLE, IERR, NAME) 

These  arguments  have  EXINTO  write  679  points  of  TIME  array  data  with  a 
title  of  LFN  TEMP03. 

EXEMPT  will  read  two  types  of  tape  format.  The  first  type  was  pre- 
viously discussed.  The  second  type,  EXEMPT  old  format,  may  be  read 
only.  The  first  record  will  contain  the  number  of  points  and  title 
information. 

An  old  format  EXEMPT  data  block  is  defined  as  follows: 

a.  Frequency  Domain  Data  (Four  Records): 

1.  Record  1-13  words  long;  (1)  one  word  integer  (N) 

indicating  number  of  words  in  each  of  the  following 
records.  (2)  12  words  of  title  information  in  CDC 

display  code. 

2.  Record  2 - N words  long  containing  the  real  part  or 
magnitude  of  a frequency  function. 

3.  Record  3 - N words  long  containing  the  imaginary  part 
or  phase  of  a frequency  function. 

4.  Record  4 - N words  long  containing  the  array  of  fre- 
quencies to  match  the  function  of  the  previous  two  records. 

b.  Time  Domain  Data  (Three  Records): 

1.  Record  1-13  words  long;  same  as  frequency  data. 

2.  Record  2 - N words  long  containing  the  values  of  a time 
function. 

3.  Record  3 - N words  long  containing  the  time  array  cor- 
responding to  the  values  in  previous  record. 

EXEMPT  will  use  the  following  method  to  determine  whether  a data 
set  is  time  or  frequency  domain  data: 


a.  A data  set  not  previously  defined  will  be  assumed  to  be  fre- 
quency domain  data.  A data  set  that  has  been  previously 
defined  will  retain  its  domain  type.  (To  classify  a data  set 
see  the  TYPE  command . ) 

b.  To  avoid  reading  wrong  data,  an  EXEMPT  error  will  occur  when 
the  number  of  points  in  the  first  record  is  the  same  as  the 
length  of  the  first  record.  Also,  when  the  number  of  points 
read  do  not  agree  with  the  number  of  points  specified,  an 
EXEMPT  error  will  occur. 

Example: 

Figure  3 is  a diagram  of  sequential  file  containing  frequency  data 
sets.  There  are  four  files  ordered  W,  Y,  X,  and  Z.  Each  data  set  on  these 
files  are  ordered  WA,  WC,  WB,  WD,  YA,  YB,  YC,  YD,  etc.  It  is  desired  to 
read  the  files  and  data  sets  alphabetically  (i.e.,  read  file  W in  the  order 
WA,  WB,  WC,  WD  followed  by  file  Y) . 

To  read  the  data  sets  in  the  order  desired,  one  must  be  able  to 
skip  files  and  blocks.  In  EXEMPT,  to  skip  n files  means  to  skip  over  n 
EOF  marks  and  position  the  head  at  the  beginning  of  the  current  file.  To 
skip  n blocks  (note:  one  block  is  one  complete  data  set)  means  to  skip 
over  n data  blocks.  Also,  skipping  zero  files  or  blocks  (i.e.,  n = 0) 
is  quite  legitimate;  however,  the  meanings  are  different.  For  blocks, 
this  is  a null  operation;  but  for  files,  the  head  is  set  at  the  beginning 
of  the  current  file. 

By  following  the  EXEMPT  command  language  flow  of  Figure  4,  one  can  see 
how  the  data  sets  are  read  in  the  desired  order. 

2.  RANDOM  ACCESS  FILES 


The  random  access  mass  storage  file  is  structured  as  a master  index/ 
subindex  file.  The  master  index  will  include  one  subindex  containing  all 
active  file  names.  There  can  be  at  most  20  user  named  files  on  the  master 
index.  For  each  file  name  there  is  a directory  containing  each  subindex 
and  its  header  whose  location  is  the  first  subindex  of  the  file.  A file 
has  54  other  sub indexes  that  are  used  as  data  array  locations.  The 
header  for  each  data  array  is  50  words  long  and  contains  the  data  array 
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REWIND  31 

READ  WA  FILE  =*  31 

SKIPB  31,  1,  FRQFIL 
READ  WB  FILE  = 31 

SKIPB  31,-2  , FRQFIL 

READ  WC  FILE  * 31 


SKIPB 

31, 

1 

, FRQFIL 

READ 

WD 

FILE  » 31 

SKIPF 

31, 

2 

READ 

XA, 

X3, 

XC,  XD 

FILE  - 

31 

SKIPF 

31, 

-1 

READ 

ya, 

Y3, 

YC,  YD 

FILE  - 

31 

SKIPF 

31. 

2 

READ 

ZA, 

ZB, 

ZC,  ZD 

FILE  = 

31 

Figure  4.  EXEMPT  Command  Language  Flow 


name,  12  words  of  title  information,  the  number  of  points,  subindex 
location  for  the  data  array  response  type  (i.e. , time,  frequency,  or  no 
type),  subindex  location  for  the  independent  data,  the  format  type,  and 
31  words  of  blank  information.  Figure  5 shows  the  file  structure  of 
the  random  access  data  base.  NOTE:  In  discussions  in  this  document, 
the  individual  subfiles  are  termed  files  and  are  referred  to  as  separate 
files.  In  fact,  they  are  subfiles  of  the  one  random  access  mass  storage 


FORTRAN  file  unit  99. 

EXEMPT  keeps  a list  of  names  associated  with  the  20  files.  This 
set  of  names  plus  the  number  of  active  files  is  written  as  the  21st 


record.  Thus,  the  21st  record  contains  only  21  words  which  are  the 
active  file  names  and  the  number  of  active  files.  When  a new  file  is 
added,  the  first  zero  entry  in  the  active  file  namelist  is  used  for  the 
index  for  the  subfile  so  that  a file  will  be  created  corresponding 
to  the  index  in  the  list  of  active  file  names. 

The  random  access  file  is  opened  with  the  master  index  in  effect. 
The  master  index  is  kept  in  the  common  block/MSTRG/  and  a subindex  cor- 
responding to  the  proper  active  file  name  is  read  from  the  random  access 
file.  Only  one  subindex  is  contained  within  EXEMPT  at  any  time.  Thus, 
when  referencing  other  files  in  the  subindex  the  file  currently  in  use 
is  rewritten  onto  the  random  access  file  before  a new  subindex  is  read. 


Particular  care  is  made  in  using  the  random  access  disk  in  this 
master  index/subindex  form  to  always  write  the  current  subindex  onto  the 
record  indicated  from  the  master  index  before  a new  subindex  is  read. 

In  this  way,  the  current  values  for  the  record  pointers  are  maintained 
correctly  so  that  data  are  not  lost  on  the  random  access  file  because  of 
an  index  error.  Each  file  is  organized  in  the  same  way,  containing  the 
file  header  and  directory  and  then  the  block  records  for  the  54  data 
arrays  on  each  file. 

The  first  subindex  is  always  used  for  the  temporary  file.  When 
opening  the  mass  storage  file,  the  temporary  file  is  always  created 
first  and  then  the  file  that  the  user  desires.  If  the  user  is  accessing 
the  temporary  file,  the  program  does  not  open  another  file.  However,  if 
the  user  has  referenced  a particular  file  by  name  or  number,  the  user's 
file  will  be  opened  as  the  second  file  using  subindex  2.  This  process 
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DATA  ARRAY 
NAME 


> 


continues  as  the  user  opens  new  files  until  all  20  files  are  in  use. 

When  the  user  requests  that  a file  be  purged,  two  procedures  are  followed. 
If  the  entire  file  is  to  be  unloaded,  the  subindex  pointing  to  the  file 
is  zeroed,  the  name  in  the  active  file  namelist  is  zeroed,  the  number  of 
files  active  is  decremented,  and  the  resulting  subindex  and  new  list  of 
file  names  are  written  onto  the  random  access  disk.  The  second  procedure 
occurs  when  the  file  is  to  be  set  empty.  In  this  case,  the  file  block 
directory  is  changed  to  indicate  that  all  the  blocks  are  empty.  Only 
the  file  header  and  directory  information  are  rewritten.  Tne  subindex, 
which  is  in  use  after  these  operations,  is  then  written  onto  the  disk 
and  operation  continues.  The  next  time  that  particular  file  is  accessed, 
the  program  recognizes  that  the  file  is  empty  and  the  first  block  will 
be  used.  Purge  file  operations  are  an  irreversible  process.  If  the 
user  indicates  a purge  file  operation,  the  contents  of  that  file  are 
removed  completely  (see  PURGE  command.  Section  II-B) . 

Lastly,  when  the  particular  EXEMPT  job  has  been  completed  and  a 
STOP  is  executed,  the  FORTRAN  I/O  operations  automatically  write  the 
master  index  to  the  end  of  the  file.  This  means  that  if  no  FORTRAN  file 
errors  occurred  during  operation,  the  random  access  file  can  be  cataloged 
as  a permanent  file  on  the  Air  Force  Weapons  Lab  computer  system.  When 
EXEMPT  completes  operation,  the  master  index  is  always  restored  and  the 
last  used  subindex  is  written  onto  the  file.  If  an  error  occurs  that  is 
not  associated  with  a FORTRAN  file  operation,  the  SCOPE  operating  system 
automatically  calls  a routine  contained  in  EXEMPT  called  CLOSRM.  No 
attempt  is  made  to  change  the  index  to  the  master  index  and  close  the 
random  access  file  correctly.  If  a fatal  error  is  detected  by  the 
FORTRAN  I/O,  the  routines  associated  with  the  FORTRAN  I/O  will  automat- 
ically try  to  close  all  FORTRAN  files.  This  means  that  the  index  in 
use  at  the  time  of  the  fatal  error  will  be  written  as  the  last  record  on 
the  random  access  file.  If  this  happens  to  be  a subindex,  then  the 
master  index  and  the  file  structure  is  completely  lost  due  to  the  error. 
This  is  a systems  function  with  no  user  control  to  override  the  FORTRAN 
error  processing.  Consequently,  it  is  not  recommended  that  the  random 
access  file  be  used  as  a permanent  file.  It  can  be  temporarily  used, 
but  a backup  must  be  maintained.  This  can  be  done  through  the  permanent 
data  master  tape. 


•j 
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The  permanent  data  master  tape  is  used  to  copy  the  random  access 
disk  to  a sequential  file  which  can  be  saved  either  as  a permanent  file 
or  as  a magnetic  tape  backup.  This  particular  file  is  written  in  such  a 
way  as  to  minimize  the  storage  requirements  for  the  file  and  yet  at  the 

same  time  allow  restoration  of  the  random  access  file  to  the  exact 

condition  at  the  time  when  the  master  tape  was  written.  The  format  for 

the  tape  structure  is  found  in  Table  4.  The  first  record  on  the  master 

data  tape  contains  the  active  file  names  and  the  number  of  active 
files.  This  is  the  same  record  written  on  the  random  access  file  as 
record  21  of  the  master  index.  This  allows  the  same  active  files  to  be 
regenerated  in  the  same  indexing  originally  obtained.  Each  active  file 
is  then  written  on  the  master  data  tape.  Each  file  is  formatted  on  the 
master  tape  by  writing  the  file  header  and  directory,  and  all  data  arrays 
referenced  in  the  directory.  Only  those  blocks  containing  data  will  be 
written  onto  the  master  tape  and  each  block  will  be  written  as  one  record. 
Frequency  dependent  data  will  contain  both  real  and  imaginary  data.  If 
the  block  is  empty,  the  data  arrays  are  not  written.  This  minimizes  the 
storage  requirements  for  the  permanent  master  tape.  In  addition,  only 
the  number  of  points  indicated  on  the  file  are  written  onto  the  master 
data  tape.  After  all  the  records  are  copied  from  one  of  the  random  access 
files,  the  next  active  file  is  found  and  the  records  for  this  file  are 
written  beginning  with  the  file  header  and  the  directory  which  contains 
2715  words,  followed  by  the  data  arrays.  It  should  be  noted  that  every 
block  header  and  directory  from  a file  is  written  on  a master  data  tape.  t 
This  allows  easy  regeneration  of  the  random  access  file  at  the  beginning 
of  execution  of  EXEMPT. 

The  master  data  tape  can  be  obtained  at  the  end  of  execution  using 
the  EXEMPT  command  SAVRAN  and,  when  it  is  written,  the  data  tape  is 
written  onto  a file  name  specified  by  the  user  or  the  default  name,  NEWDAT. 
At  the  beginning  of  execution,  a previously  existing  master  data  tape  can 
be  read  into  the  random  access  file  at  the  time  that  the  random  access 
file  is  opened  by  the  EXEMPT  command  GETRAN  and  requesting  the  user's  named 
file  or  the  default  name  OLDDAT  for  the  master  data  tape.  The  master 
data  tape  is  automatically  unloaded  by  the  EXEMPT  program  after  the  random 
access  file  is  generated. 
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RECORD 

1 

2 

3 

4 

5 


Table  4 

PERMANENT  DATA  MASTER  TAPE  STRUCTURE 

NO.  OF 

WORDS  CONTENTS 

21  ACTIVE  FILE  NAMES,  NUMBER  OF  ACTIVE 

FILES 

2715  FILE  HEADER  AND  DIRECTORY  FOR 

TEMPORARY  FILE 

N DATA  ARRAY  1 FOR  TEMPORARY  FILE 

N DATA  ARRAY  2 FOR  TEMPORARY  FILE 

N DATA  ARRAY  3 FOR  TEMPORARY  FILE 


J 

J+l 

J+2 


2715  FILE  HEADER  AND  DIRECTORY  FOR 

NEXT  ACTIVE  FILE 

M DATA  ARRAY  1 FOR  NEXT  ACTIVE  FILE 

M DATA  ARRAY  2 FOR  NEXT  ACTIVE  FILE 

FOR  BLOCK  1 
DATA  FOR  BLOCK  1 
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The  random  access  data  base  will  reside  on  LFN  TAPE99.  If  the  user 
has  an  existing  EXEMPT  compatible  random  access  data  base  he  wishes  to 
use,  it  must  be  attached  with  a LFN  of  TAPE99.  Likewise,  to  save  an 
EXEMPT  random  access  data  base  for  later  use,  a LFN  of  TAPE99  must  be 
used.  However,  if  the  user  has  attached  a random  access  data  base  on 
LFN  TAPE99,  he  must  use  the  EXTEND  or  ALTER  option  at  the  end  of  the 
EXEMPT  run  in  order  to  save  the  modified  data  base.  As  mentioned  previously, 
it  is  recommended  that  the  permanent  master  data  tape  system  be  used. 

EXEMPT  has  the  capability  of  writing  binary  data  files  which  may  be 
saved  as  external  data  tapes  if  required.  In  addition,  specific  files, 
in  particular  the  master  data  discussed  earlier  in  this  subsection,  are 
written.  Any  of  the  sequential  data  capes  written  using  the  standard 
EXEMPT  sequential  file  commands  can  be  requested  as  an  external  tape  and 
saved.  The  file  names  associated  with  these  tapes  are  the  names  given 
by  the  user  unless  they  are  specified  with  a number  between  1 and  99.  In 
that  case,  the  file  name  will  be  TEMPXX  where  the  number  is  coded  into 
the  position  XX.  The  tens  digit  is  not  removed  for  the  numbers  between 
1 and  9 so  that  these  files  become  TEMP01  through  TEMP09.  The  random 
access  data  file  is  saved  on  disk  and  can  be  cataloged  as  a permanent 
file,  but  cannot  be  used  with  a magnetic  tape  copy.  In  order  to  save  this 
information  permanently,  the  user  should  use  the  GETRAN/SAVRAN  option 
to  write  a permanent  data  master  tape. 
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SECTION  IV 

USER  INTERFACE  REQUIREMENTS 


The  user  may  interface  data  directly  with  EXEM  PT  by  using  the 
rigid  interface  specification  provided  by  the  OPRn,  SUBa,  MYCOMn,  and 
NAME  CALLS  discussed  in  the  ECL  of  Section  II.  In  addition,  data  from 
external  or  internal  programs  may  be  retrieved  from  or  written  to  a 
peripheral  sequential  file  using  the  format  specified  in  Section  III. 

When  using  the  OPRn,  SUBa,  MYCOMn,  or  NAME  subroutine  call,  the 
interface  must  be  as  specified  below: 

Subroutine  Interface 


) 


Name 

SUBROUTINE  ) OPRn  ( (RSP,M,N,X,C[  ,A] , B]  ]) 

SUBa  l 
MYCOMn  / 

RSP  - Global  real  array  of  up  to  20  values  used  as  parameters 
in  the  subprogram 

M Maximum  number  of  points  found  in  X,  C,  A,  or  B 

N Specifies  type  of  resultant 

1 Time  or  no  type 

2 Frequency 

X Independent  array 

C Resultant  array 

A Optional  input  array 

B Optional  input  array 


DIMENSION 

RSP(20) 

LEVEL  2,  X,  C [,  A [, 

B]] 

DIMENSION 

X(M) 

DIMENSION 

C (M,N) 
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All  arrays,  except  RSP,  are  passed  via  the  array  BLNKCM  located  in 
blank  common.  BLNKCM  is  declared  a LEVEL  2 variable  so  that  it  mav  be 
loaded  in  LCM  on  the  CDC  7600.  As  such,  all  other  variables  that 
reference  any  address  in  the  array  BLNKCM  must  also  be  declared  LEVEL  2. 
On  the  CDC  6600,  the  LEVEL  2 declaration  is  ignored. 

When  complex  data  is  passed  through  the  user  subroutine  or  to  be 
returned  through  the  array  C,  the  first  column  (i.e.,  C(M,1),  A(M,1),  and 
B(M, 1) ) will  contain  the  real  data  and  the  second  column  (i.e.,  C(M,2), 
A(M, 2) , and  B(M,2))  will  contain  the  imaginary  data. 

The  user  must  replace  the  EXEMPT  dummy  subroutines  OPRn,  SUBa,  and 
MYCOMn  on  the  relocatable  binary  file  by  using  the  COPYL  SCOPE  utility. 

To  use  the  NAME  subroutine  call,  the  user  must  append  the  entry  point 
identified  as  NAME  to  the  Direct  Manipulation  Processor  segment  using 
the  COPYL  with  the  append  option  specified. 

To  append  the  entry  point  for  a NAME  subroutine,  the  user  must 
replace  the  EXEMPT  subroutine  CNGSUB  on  the  relocatable  binary  file  by 
using  the  COPYL  SCOPE  utility.  The  replacement  must  include  an  EXTERNAL 
declaration  followed  by  a call  to  the  EXEMPT  subroutine  SETUP.  This 
sequence  will  add  the  user's  subroutine  name  and  address  location  to  a 
table  contained  within  EXEMPT.  The  arguments  to  SETUP,  as  shown  below. 


include  the  user  subroutine  name  in  Hollerith  form  followed  bv  the 


subroutine  name.  As  many  pairs  of  arguments  are  required  as  there  are 
user-named  subroutines,  up  to  a maximum  of  36  pairs. 


SUBROUTINE  CNGSUB 

EXTERNAL  PULSE1 

CALL  SETUP (6HPULSE1, PULSE1) 

RETURN 

END 


1. 


NO  INPUT  ARGUMENTS  OR  RESULTANT  SPECIFIED 


| 

The  MYCOMn  subroutines  should  be  che  only  subroutines  callable  In 
this  form.  However,  since  the  interface  is  the  same  for  all  user  subroutine 
calla,  any  subroutine  type  may  be  called  in  this  form.  The  MYCOMn 
subroutines  will  be  the  only  user  subroutines  that  may  have  an  alternate 
name  assigned.  For  example,  consider  Che  following  ECL: 

MYC0M1  - GTD 
CALL  GTD(  ) 

The  MYCOMn  subroutines  are  Intended  to  be  used  as  alternate  user- 
defined  commands.  Since  this  is  a transfer  of  control,  che  CALL  keyword 
must  be  used  in  conjunction  with  the  alternate  command.  Also,  note  the 
use  of  the  null  parentheses.  These  must  be  present  since  they  delimit 
the  start  and  end  of  a subroutine  parameter  list,  even  if  a null  list. 

The  interface  for  the  above  example  and  all  null  parameter  list 
calls  is: 

SUBROUTINE  MYC0M1 (GREAL,  M,  N,  X,  C) 
where: 

GREAL  - Global  real  array 

M - Set  to  zero 

N - Set  to  zero 

X - Dummy  argument 

C - Dummy  argument 

2.  RESULTANT  ONLY  SPECIFIED 

All  subroutines  may  be  used  in  this  manner.  However,  there  are 
three  types  of  interfaces  within  this  class  and  all  have  to  do  with  how 
the  independent  is  to  be  passed. 


u 


a.  If  the  resultant  is  a frequency  or  time  domain  response,  the 
resultant’s  associated  independent  is  called  through  the 
subroutine  interface. 

b.  If  the  resultant  is  classified  as  time  or  frequency  (e.g., 

TYPE  command)  and  the  resultant  has  not  yet  been  defined,  the 
currently  defined  frequency  or  time  independent  is  called 
through  the  subroutine  interface. 

c.  If  the  resultant  is  not  classified  as  frequency  or  time,  a 
null  array  is  passed  as  the  independent  with  the  following 
input  parameters: 

1.  M - Set  to  1/3  of  current  BLNKCM  size 

2.  N - Set  to  2.  On  return,  user  sets  N to  the  following: 

N ■ 0,  resultant  is  a no  type 

N » 1,  resultant  is  a time  function 
N - 2,  resultant  is  a frequency  function 

3.  X - On  return  contains  independent  values  for  N f 0 

The  following  are  examples  of  calls  to  subroutines  without  optional 
input  arrays. 

Example  1: 

A - SUBA(  ) 

A is  the  resultant  data  set  for  the  subroutine  SUBA.  A will  be 
classified  to  be  frequency  or  time  domain  data  or  a no  type  data  by  the 
user  when  the  SUBA  subroutine  call  is  completed. 

Example  2: 

TIME  - LIN101 
B - 1. /TIME 
LIST  B 

TIME  - LIN51 
CALL  0PR1(B) 

B is  the  resultant  data  set  for  the  subroutine  0PR1.  Since  B had 
been  classified  as  a time  domain  response  in  statement  2,  the  independent 
associated  with  B,  statement  1,  is  called  through  the  interface  to  0PR1. 
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3.  ONE  OR  TWO  INPUT  ARRAY  AND  RESULTANT  SPECIFIED 

There  are  three  types  of  subroutine  Interfaces  within  this  class 
and,  again,  all  have  to  do  with  the  way  the  independent  is  to  be  passed: 

a.  Same  as  the  no  input  array  interface  (1).  Unless  TIME  or  FREQ 
is  specified  as  the  first  input  array,  then  the  specified 
global  independent  is  passed  through  the  interface  and  will 
become  the  independent  data  associated  with  the  resultant. 

b.  Same  as  the  no  input  array  interface  (2).  Unless  TIME  or  FREQ 
is  specified  as  the  first  input  array,  then  the  specified 
global  independent  is  passed  through  the  interface  and  will 
become  the  independent  data  associated  with  the  resultant. 

c.  If  the  resultant  is  not  typed,  the  independent  associated 
with  the  input  array  is  called  through  the  interface. 

d.  If  the  input  array  is  not  typed  then  a dummy  independent  of 

the  same  size  of  the  input  array  is  passed.  Also,  the  resultant's 
array  is  set  to  the  size  of  the  input  array.  However,  since 
the  interface  does  not  know  at  this  point  whether  the  resultant 
will  be  frequency  or  time  domain,  a doubly  dimensioned  resultant 
array  will  be  passed  (i.e.,  dimension  C to  (M,  2)).  On  the 
return,  M will  be  the  number  of  points  in  independent  and/or 
resultant.  N will  be  the  type  of  resultant  generated: 

1.  N - 0 - No  Independent  Associated 

2.  N - 1 - Time  Independent 

3.  N ■ 2 - Frequency  Independent 

4.  PASSING  THE  GLOBAL  REAL  ARRAY  AND  CURRENT  INDEPENDENT  ARRAY 


I 


The  global  real  array  is  set  externally  to  the  calling  sequence 
through  the  variables  R1-R20.  These  variables  are  set  prior  to  the  sub- 
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TIME  = LIN128 , 0.5E-6 

R2  = 0 $ RESET  GREAL(2)  = 0. 

R3  ■ 9.366E6  $ SET  GREAL  (3)  - 9.366E6 

R4  - 1.953E7  $ SET  GREAL  (4)  = 1.953E7 

TYPE  Y » TIME  $ CLASSIFY  Y AS  A TIME  FUNCTION 

CALL  PULSEl(Y) 

LIST  Y 
END 

In  the  above  example,  the  user  may  wish  to  pass  the  TIME  array  and 
set  Y internally  as  a time  function.  In  this  case,  the  TYPE  command 
would  be  omitted  and  the  subroutine  call  would  be: 

CALL  PULSE1 (Y, TIME) 

Also,  on  return  from  subroutine  PULSE1  the  user  need  not  set  N * 1 
to  signify  a time  function.  EXEMPT  will  do  this  on  seeing  TIME  in  the 
calling  sequence.  Additionally,  the  call  could  have  been  placed  in  an 
arithmetic  statement  such  as: 

Y « PULSEl(TIME) 

OR 

TYPE  Y = TIME 

Y = PULSE1(  ) 

A complete  example  of  passing  the  global  real  array  and  the  subroutine 
interface  is  shown  in  Figures  6 and  7. 
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THIS  PAGE  IS  BEST  QUALITY  FRACTICABLI 
from  OQPY  PUKMISHEL  TO  HDC 


SUHROUTINE  PUL SE 1 ( RSP , H , N , T » r T ) 
DIMENSION  PSP(IO) ,T(1) ,FT(1) 

A = RSP(2) 

ALPHA  = RSp ( 3 ) 

F = RSP ( A ) 

Q = RSP ( 5 ) 

TSHIFT  = RSP ( 6 ) 

NT  = M 

DATA  TW0PI/6. 23318631/ 

W ~ T wOP I *F 
IF (O.EO.O.O)  30  TO  1 
ALPHA  = 0 . 5*w / 3 

1 CONTINUE 

IF(A.NE.O.O)  30  TO  2 
TP  = ATAN (W/ALPHA ) /* 

A = 1 . 0/ (EXP (-ALPHA*TP) -SIM (W*TP)  ) 

2 CONTINUE 

DO  10  I = 1* NT 
X = T ( I ) ♦ TSHIFT 

r T ( I ) = A^EXP  ( -ALPHA*X  ) *SIN  (V/*X  ) 

10  CONTINUE 

X = - TSHIFT 
I = 0 

15  I = I ♦ 1 

IF (T ( I ) . GE . X ) GO  TO  20 
FT ( I ) = 0.0 
GO  TO  15 
20  CONTINUE 
RETURN 
END 


Figure  7.  User  Subroutine  Interface 
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SECTION  V 

COMPUTER  INTERFACE  REQUIREMENTS 


There  are  three  permanent  filer  with  which  the  user  will  come  in  contact. 

a.  EXEMPT  Binary  File 

b.  EXEMPT  Segmentation  Directive  File 

c.  EXEMPT  Absolute  File 

These  files  are  more  clearly  defined  in  the  EXEMPT  Programmer's  Manual.  At 
this  point,  the  user  only  needs  to  know  that  the  binary  and  segmentation 
directive  files  will  be  needed  when  using  the  COPYL  utility.  The  user 
should  normally  come  into  contact  with  only  the  absolute  file.  The  following 
examples  should  illustrate  the  point. 

1.  EXAMPLE  1 DESCRIPTION 

Example  1,  shown  in  Figure  8,  illustrates  some  of  the  EXEMPT 
arithmetic  capabilities  available.  The  example  computes  the  transfer 
function 


H(u>) 


(a  + jw) 

for  each  frequency  value.  The  output  listing  is  shown  in  Figure  9. 


EXAMPLE  2 DESCRIPTION 


Example  2 calculates  the  input  impedance  of  a shorted  transmission 
line  by  using  some  of  the  standard  functions  available  within  EXEMPT.  The 
ECL  is  shown  in  Figure  10,  the  LIST  output  is  shown  in  Figure  11,  and 
the  resultant  METAPLOT  is  shown  in  Figure  12. 

For  a shorted  transmission  line,  the  input  impedance  is  given  by 

Z__  - Z tanh  (y2.) 

SC  o 

However,  since  hyperbolic  functions  are  not  included  among  the 
standard  functions  available  to  EXEMPT,  the  exponential  function 


S3 


i 


J08  CARO. 

ACCOUNT  CARO. 

ATTACH. ABS.PFNAME  (PARAMETER  LIST). 

ABS. 

Si  7/8/9  (EOR)  EXEMPT  ECL 

FREQ  = LIN51.50.0E6 
ALPHA  = 10.0E6 
JH  * CMPLX ( 0 . 6.2831 8*FREQ) 

HH  r 1/ (ALPHA  ♦ JW) 

LIST  HW 
ENO 

* 6/7/B/9  (EOF)  ENO  OF  JOB 


EXEMPT  ABSOLUTE  FILE 


Figure  8.  Job  Stream  and  ECL  for  Example  1 
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LIST  OF  DATA  SET  HH 


INOEX 

FREQ 

<466 

PHASE 

1 

0. 

1.000E-07 

0. 

2 

1 ♦OOOE*06 

8.467E-08 

-32.1 

3 

2.000E*06 

6.227E-08 

-51.5 

4 

3.000E+05 

4.687E-08 

-62.1 

5 

4. 000E+06 

3.697E-O0 

-68.3 

6 

5.000E*05 

3.033E-08 

-72.3 

7 

6.000E+06 

2.564E-O0 

-75.1 

8 

7.000E*06 

2.217E-08 

-77.2 

9 

6. 000E  + 0S 

1.951E-08 

-78.7 

10 

9.000E*06 

1.741E-08 

-80.0 

11 

1 • 000E*07 

1.572E-08 

-81.0 

12 

1 • 1 00E*07 

1.432E-08 

-81.9 

13 

1.200E*07 

1.315E-08 

-82.4 

14 

1.300E*07 

1.215E-Q8 

-83.0 

15 

1.400E*07 

1.130E-03 

-83.5 

16 

1.500E*07 

1.055E-08 

-83.9 

17 

1.600E+07 

9.898E-09 

-84.3 

18 

1.700E*07 

9.321E-09 

—84 . 7 

19 

1.800E*07 

8.808E-09 

-84.9 

20 

i.900E*07 

8.347E-09 

-85.2 

21 

2.000E«-07 

7 « 933E-09 

-85.5 

22 

2.100E*07 

7.557E-09 

-85,7 

23 

2.200E*07 

7.21SE-09 

-85.9 

24 

2.30QE*07 

6.903E-Q9 

-86.0 

25 

2.400E*07 

6.617E-09 

-86.2 

26 

2«500E*07 

6.353E-09 

-86.4 

IMDEX 

FREQ 

HAG 

PHASE 

27 

2.600E*07 

6.110E-09 

-86.5 

26 

2. 700E+07 

5.884E-09 

-86.6 

29 

2.800E+07 

5.675E-09 

-66.7 

30 

2.900E*07 

5.480E-09 

-86.9 

31 

3.000E*07 

5.298E-09 

-87.0 

32 

3.100E*07 

5.127E-09 

-87.1 

33 

3.200E*07 

4.967E-09 

-87.2 

34 

3.300E+07 

4.817E-09 

-67.2 

35 

3.400E*07 

4.676E-09 

-87.3 

36 

3.500E+07 

4.543E-09 

-87.4 

37 

3.600E*07 

4.417E-C9 

-87.5 

38 

3.700E*07 

4.299E-09 

-87.5 

39 

3.800E*07 

4.18SE-09 

-87.6 

40 

3.900E*07 

4.078E-09 

-87.7 

41 

. 4.000E*07 

3.976E-09 

-87.7 

42 

4.100E*07 

3.879E-09 

-87.8 

43 

4.200E*07 

3.787E-09 

-67.8 

44 

4.300E*07 

3.699E-09 

-87.9 

45 

4.400£*07 

3.615E-09 

-87.9 

46 

4.500E*07 

3.535E-09 

-68  . 0 

47 

4.600E*07 

3.458E-09 

-68.0 

48 

4. 700E*07 

3.384E-09 

-68.1 

49 

4.800E*07 

3.314E-09 

-68.1 

50 

4.900E*07 

3.246E-09 

-68.1 

51 

5.000E«-07 

3.181E-09 

-88.2 

ENO  OF  LIST  • 


Figure  9.  Output  for  Example  1 
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JOB  CARO. 

ACCOUNT  CARO. 

REQUEST(PLTFIL**Q) 

OISPOSE(PLTFIL»*MF*PCF)  dispose  plots  to  calcomp 

ATTACH* A8S*PFNAM£  (PARAMETER  LIST).  EXEMPT  ABSOLUTE  FILE 

ABS. 

& 7/8/9  (EOR)  EXEMPT  ECL 

S 

S IMPEDANCE  OF  SHORTED  TRANSMISSION  LINE 

S 

PI  * 3.141592654 
FREQ  s 1.E6*LIN50*50.E6 
« a 2*PI*FREQ 


LENGTH  a 2 
ZO  a 50 
EPSR  a 3 
ROC  a 0.1E-3 
C a 3.0E8 


S CHARACTERISTIC  IMPEDANCE 
S RELATIVE  DIELECTRIC  CONSTANT 
S DC  RESISTANCE  PER  METER 
S SPEED  OF  LIGHT 


GL  a CMPLX{RDC*SQRT(-REQ)/(2*Z0) *W*SQRT (EPSR) /C) *LENGTH 
S GAMMA  X LENGTH 

ZSC  a Z0*((EXP(GL)  - EXP(-GL) )/(EXP(GL)  ♦ EXP(-GL))) 

S ZSC  * ZO*TANH(GAMMA*LENGTH) 

S TANHi(X)  a ( (EXP  (X)  -EXP (-X)  ) / (EXP  ( X)  *EXP (-X)  ) ) 

PLOT  ZSC  TITLEY  a *IMPEDANCE* 

list  GL.ZSC 

END 

• 6Z7/8/9  (EOF)  END  OF  JOB 


Figure  10.  Job  Stream  and  ECL  for  Example  2 
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CIST  OF 

DATA  SET 

OL 

INOEX 

FREU 

RAG 

PHASE 

INDEX 

FREQ 

MAG 

Phase 

1 

1 < OOOE  *06 

7.258E-02 

88.4 

26 

2.600E*07 

1.89 

89.7 

2 

2.000E*06 

.145 

88.9 

27 

2. 700E*07 

1.96 

89.7 

3 

3.000E*06 

.218 

89.1 

28 

2.800E*07 

2.03 

89.7 

* 

4. 000E  *06 

.290 

89.2 

29 

2.900E*07 

2.10 

89.7 

S 

5.000E*06 

.363 

89.3 

30 

3.000E*07 

2.18 

89.7 

6 

6.000E*06 

.435 

89.4 

31 

3. 100E*07 

2.25 

89.7 

7 

7.000E*06 

.508 

89.4 

32 

3«200E*07 

2.32 

89.7 

8 

8.000E«06 

.580 

89.4 

33 

3.300E*07 

2.39 

89.7 

9 

9.000E*06 

.653 

89.5 

34 

3.400E*07 

2.47 

89.7 

10 

1 •OOOE *07 

.726 

89.5 

35 

3.500E*07 

2.54 

89.7 

11 

1.100E*07 

.798 

89.5 

36 

3.600E*07 

2.61 

89.7 

12 

1.200E*07 

.871 

89.5 

37 

3.700E*07 

2.68 

89.7 

13 

1.300E*07 

.943 

89.6 

38 

3.800E*07 

2.76 

89.7 

1* 

1 • 400E  *07 

1.02 

89.6 

39 

3.900E*07 

2.83 

89.7 

15 

1 .500E  *0  7 

1.09 

89.6 

40 

4.000E*07 

2.90 

89.8 

16 

1.600E*07 

1.16 

89.6 

41 

4. 1 OOE*  07 

2.97 

89.8 

17 

1.700E»07 

1.23 

89.6 

42 

4.200E*07 

3.05 

89.8 

18 

1.800E*07 

1.31 

89.6 

43 

4.300E*07 

3.12 

89.8 

19 

1 • 900E»0  7 

1.38 

89.6 

44 

4.400E*07 

3.19 

89.8 

20 

2.000E*07 

1.45 

89.6 

45 

4.500E*07 

3.26 

89.8 

21 

2.100E*07 

1.52 

89.7 

46 

4.600E*07 

3.34 

89.8 

22 

2.200E*07 

1.60 

89.7 

47 

4. 700E*0  7 

3.41 

89.8 

23 

2.300E*07 

1.67 

89.7 

48 

4.800E*07 

3.48 

89.8 

2* 

2.400E*07 

1.74 

89.7 

49 

4.900E*07 

3.56 

89.8 

25 

2.500E*07 

1.81 

89.7 

50 

5.000E*07 

3.63 

89.8 

LIST  OF 

DATA  SET 

ZSC 

INOEX 

FREQ 

RAG 

PHASE 

INOEX 

FRE  a 

MAG 

PHASE 

1 

1. OOOE *06 

3.64 

88.4 

26 

2.600E*07 

153. 

-88.0 

2 

2.000E*06 

7.31 

88.9 

27 

2.700E.07 

122. 

-88.3 

3 

3.000E*06 

n.i 

89.1 

28 

2.800E*07 

101. 

-88.5 

6 

4.00 OE *06 

14.9 

89.2 

29 

2.900E*07 

84.7 

-88.6 

5 

5 . OOOE • 06 

19.0 

89.2 

30 

3.000E*07 

72.2 

-88.7 

6 

6. OOOE *06 

23.3 

89.3 

31 

3. 1 OOE  *07 

62.0 

-88.7 

7 

7. OOOE *06 

27.8 

89.3 

32 

3.200E*07 

53.6 

-88.7 

8 

8. 000E*06 

32.8 

89.3 

33 

3.300E*07 

46.3 

-88.7 

9 

9.000E*06 

38.2 

89.3 

34 

3.400E*07 

40.0 

-88.6 

10 

1 • OOOE  *0  7 

44.3 

89.3 

35 

3.500E*07 

34.4 

-88.5 

11 

1 . 1 OOE  *0  7 

51.3 

89.2 

36 

3.600E*07 

29.3 

-88.4 

12 

1.200E*07 

59.3 

89.2 

37 

3.700E*07 

24.6 

-88.2 

13 

1 . 300£»0  7 

58.9 

89.1 

38 

3.800E*07 

20.2 

-88.0 

14 

1.400E*07 

90.6 

89.0 

39 

3.900E«07 

16.1 

-87.6 

15 

1.500E*07 

95.4 

88.9 

40 

4.000E«07 

12.2 

-86.9 

16 

1.600E*07 

115. 

88.7 

41 

4.1 OOE*07 

8.45 

-85.5 

17 

1.700E*07 

142. 

88.5 

42 

4,200E*07 

4.78 

-82.1 

18 

1 . 800E*07 

184. 

88.1 

43 

4.300E»07 

1.27 

-59.0 

19 

1.900E*07 

257. 

87.3 

44 

4 , 400E*07 

2.62 

75.3 

20 

2.000E*07 

414. 

85.7 

45 

4.500E*07 

6.23 

83.7 

21 

2.100E*07 

1.039E*03 

79.0 

46 

4.600E«07 

9.94 

85.9 

22 

2.200E*07 

1 . 850E*03 

-69.7 

47 

4.700E*07 

13.8 

86.9 

23 

2.300E*07 

507. 

-84.4 

48 

4.800E*07 

17.8 

87.5 

24 

2.400E*07 

290. 

-86.6 

49 

4 . 900E*07 

21.9 

87.8 

25 

2.500E ♦ 07 

202. 

-87.5 

50 

5.000E*07 

26.4 

80.0 

••••  ENO 

OF  CIST 

Figure  11.  LIST  Output  for  Example  2 
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representation  of  the  hyperbolic  tangent  is  used.  Note  that  the  exponential 
function,  EXP,  can  handle  a real  or  complex  argument.  That  is,  if 


then 


but  if 


X = real  number 
EXP(X)  - eX 


then 


X - a + jfi 

EXP  (X)  =■  ea(cos6  + jsinS) 


3.  EXAMPLE  3 DESCRIPTION 


Example  3,  reference  Figure  13  for  ECL,  is  an  exercise  in  reading 
a data  tape  generated  by  another  program  and  storing  the  generated  EXEMPT 
random  file  data  base  on  permanent  files.  The  LIST  CURFIL  output,  shown 
in  Figure  14,  lists  each  data  set  as  they  occur  on  the  random  file 
including  each  independent.  In  this  case  of  the  six  data  sets,  GTDF1- 
GTDF6,  there  are  five  independents  also  listed.  There  are  only  five 
independents  listed  because  data  sets  GTDF3  and  GTDF6  reference  the  same 
independent  and,  therefore,  this  independent  is  listed  but  once.  The 
output  for  the  SAVRAN  command  is  shown  in  Figure  15.  Figure  15  is  an 
index  to  the  data  base  saved  as  a result  of  the  SAVRAN  command. 


89 


>6  W 


AD-A066  806  BOM  CORP  ALBUQUERQUE  N MEX 
EXEMPT  USERS  MANUAL. (U) 

JAN  79 


UNCLASSIFIED 

BDM/A-76-1 

40-TR-l 

12 

2 of  2 
^066806 

■ 

njiui 

tii  .ui 

!»•«* 
iii  :mi 

iii  iii 

it!.!".* 

iii  Ui 

iiTui 

IBIS  PAGE  IS  BEST  QUALITY  PRACTICABLE 
FHOM  C VrX  JfutuUamD  TO  MG  ^ 


JOB  CARO. 

ACCOUNT  CARO. 

ATTACH*  TEMPO  1 .PFNAME.  (PARAMETER  LIST).  EXTERNAL  FILE 
REQUEST .NEWDAT  * *PF • 

REQUEST »TAPE99»*pF  • 

ATTACH. ABS.PFNAME  (PARAMETER  LIST).  EXEMPT  ABSOLUTE  FILE 

ABS. 


CATALOG. TAPE99.PFNAMEI  (PARAMETER  LIST).  DATA  BASE  SAVED  ON  RANDOM  FILE 

CATALOG.NEWDAT.PFNAMEi  (PARAMETER  LIST).  DATA  BASE  SAVEO  ON  SEQUENTIAL  FILE 

L 7/6/9  (EORJ  EXEMPT  ECL 

S EXAMPLE  TO  DEMOSTRATE  EXEMPT  LINKAGE  TO  A FILE  GENERATED  8V 
S ANOTHER  PROGRAM*  IN  EXEMPT  FORMAT. 

S ALSO*  THE  GENERATED  DATA  BASE  IS  SAVEO  ON  A PERMANENT  FILE  FOR 
S LATER  USE. 

REWIND  1 
FILE  GTDR1 
READ  GTDF1  FILE«1 
RE AO  GTDF2  FIlE*1 
READ  GT0F3  FIlE*1 
READ  GTDFA  FILE*! 

READ  GTDF5  FIlE*1 
READ  GTDF6  FIlE*1 
LIST  CURFIL 
SAVRAN 
END 

• A/7/8/9  (EOF)  ENO  OF  JOB 


Figure  13.  Job  Stream  and  ECL  for  Example  3 
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••••  c u 1 1 i n i rite 


I M 0 I * 


•••• 


THIS  PAGE  IS  BEST  QUALITY  FRACTICABIil 
miM  OOi  I ITUKWiSiliiU  TO  UDC  - ' 


rue  mocx*  t rue  nanc*  *tosi 
data  mt  nut 


mi  «TDM  run-ami* 

DATA  SIT  NANC 

•ton 


• T4CTA*  ••  OOSCRVATION  AOINT-I*  IIS.  ANAL C A I.  Jl  AMI  MOOt 
RCSAONSC  TVAC  NUN  Of  ATS  INOCX  I NO  IHOtl 

arcqucncy  Ik  l l 


00.16. AT.  »l/tS/TT  DATA  A 


LIST  or  DATA  SCT 


OTOfl 


•TOM  RUN-AMI*  • THCTA*  ••  OBSCRVXTION  AOINT-I*  Its.  AN«te*  0.  Jl  AMI  NOOC  ••.!«. AT.  »1/»S/TT  0-1  DATA  A 


INOCX 

A SCO 

RIAL 

|NAO  INOCX  ARCO 

ACAL 

IMA# 

1 

-J.kTIC-OS  ••.IS«C-«S 

SO.O 

IA  -I.6AOC-OA 

•3 

S30C-0S 

1(0. 

1 

6.00IE-OS 

k.031£*0S 

AO.O 

IS  t.ll«C-0S 

mo. 

3 

-••kOTC-OS  -S.kAIC-Ok 

70.0 

Ik  -I.10IC-0S 

ATtC-OS 

too. 

A 

j.stkc-os  -i.sAie-os 

■0.0 

it  i.oose-os 

- 

A33C-Ok 

no. 

S 

-A.IS«C-«S 

s. one-os 

•0.0 

10  -I.S1AC-0S 

•TOC -OS 

110. 

k 

-i.tioe-os  -i.iose-os 

100. 

If  k.kSAC-Ok 

AAAC-OS 

130. 

T 

I. ITSC-8S 

i.issc-ts 

110. 

10  1. SIOC-OS 

abac -os 

IAO. 

■ 

-I.AAOt-OS  -T. JA7C-0A 

110. 

It  -A.07SC-00 

■A3C-0k 

ISO. 

• 

-A.ISTC-BT  -I.AIlftS 

no. 

II  A.OOSC-Ok 

m 

OIAC-OS 

IkO. 

IS 

-l.DIC-tk 

k.Sk«C-«« 

IA0. 

i)  i.jssc-oo 

••AC-Ok 

ITO. 

II 

-).tA*C-SS  -T.AOSC-OT 

ISO. 

IA  -I.OOIC-OS 

• 

ooie-OA 

1(0. 

II 

I.StAC-SS  -1.0B3C-0S 

IkO. 

IS  I.SA3C-0S 

• 

|AAC-*T 

1*0. 

1) 

-S.SSSC-Ok 

2.AAIC-0S 

ITO. 

to  -A.askC-oo 

1 

AAIC-OS 

300. 

DATA  SCT  TITLC 


DATA  SCT  NAMC 


RCSAONSC  TV At 


HUM  or  ATS 


INDIA 


INO  I HOC A 


ARta 

•ONC 

Ik 

1 

0 

LIST 

or  OATA  SCT 

MCO 

•OCX 

VALUC  MOCX 

VALUC 

NOCX 

VALUC 

•PCX 

VALUC 

NOCX 

VALUC 

1 

-3.A7IC-05 

1 

k.OOIC-OS 

3 

-•.kOTC-BS 

A 

3.IAAC-0S 

s 

-A.ISOC-OS 

k 

-I.010C-0S 

1 

I.1TSC-0S 

• 

-i.om-os 

• 

-A.ISTC-OT 

10 

-1.I3IC-0A 

11 

-3.0ASC-OS 

It 

I.OOAC-OS 

13 

-S.SBSC-Ok 

IA 

-l.kkOC-Ok 

IS 

t.IIAC-OS 

Ik 

-3. 10IC-0S 

IT 

t.oosc-os 

IS 

-I.SIAC-OS 

It 

k.kSAC-Ok 

10 

1.31*1-08 

It 

-A.07sc-0k 

II 

A.OkSC-Ok 

13 

1.3S0C-0B 

IA 

•1 .otsc-os 

IS 

I.SA3C-8S 

Ik 

-*.00k£-0k 

OATA  SCT 

T1UI 

tin 

•TOM  AUM-AMl* 

0 THCTA*  ••  OBSERVATION  AOINT-l*  Its.  ANOLlA  0 

. JA  AMI 

OATA  SIT  NANI 
6T0AI 


RCSAONSC  TVAC 
ARCOUCNCV 


NUN  or  ATS 
Ik 


INOCX 

s 


mo  inoca 

A 


•k.li.AT.  01/0%/TT  O-l  DATA  A 


LIST  or  OATA  SCT 


OTDft 


• TON  RUN-AMt*  • THCTA*  «t  OBSERVATION  AOINT-I*  IIS.  ANSlC*  ».  JA  AMI  MOOC  ••.!*. AT.  BI/BS/TT  0-1  DATA  A 


INOCX 

ARCO 

rial 

IMA* 

INOCX 

RCAL 

1 MAO 

1 

-5.016C-03 

•T 

.ISIC-OS 

SO.O 

IA 

-5.0IOC-03 

s 

0S7C-0A 

ISO. 

I 

-S.SA1C-03 

.IAAC-OA 

AO.O 

• | ■ 

-S.IIS£-03 

• 

10AC-0A 

mo. 

3 

-S.OI1C-OI 

.TOOC-OA 

70.0 

16 

3 

TOOL -Of 

200. 

♦ 

-A.6I9C-03 

-1 

.imc-OA 

00.0 

IT 

-4.620C-03 

Tooe-OA 

no. 

s 

-S.03VC-03 

•6 
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atON  aun-»mi»  a thcta.  «a  osscavatiom  »oint*c«  its.  anolc*  t.  « am  mooc  bb.ia.at.  bi/bb/tt  a-i 


INOCA 

rato 

Hal 

INAO  INOCA 

rata 

atAL 

IN  AO 

1 

-S. 2001*13 

0761-12 
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90.0  la 

1.67TC 

13  9.ITSC-13 

220. 

1.STSC*I3 

A03I-I2 

100.  IB 

I.530C 

ii  T.asac-ii 

230. 

T 

1.1901*13 

7301-13 

110.  20 

i.osai 
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21 

-1.0571*1 
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S.299C-03  -1.3J2C-93 
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GTDFA  FREQUENCY  26 

OATA  SET  TITLE 
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Figure  15.  Index  From  SAVRAN  Command 
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A.  EXAMPLE  4 DESCRIPTION 


Figure  16  is  the  job  stream  and  ECL  for  linking  a user-defined 
subroutine  to  EXEMPT.  If  the  user  wishes  to  use  a SUBa,  OPRn,  or  a MYCOMn 
then  the  job  stream  is  exactly  the  same.  The  SCOPE  COPYL  utility  modifies 
or  appends  any  new  routines  to  the  present  program.  The  output  of  this 
example  is  presented  in  Figures  17  to  19. 
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ACCOUNT  CARO. 

ftn. 
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A • RSP<2) 

ALPHA  ■ RSPC3I 
F ■ RSP < 4 1 
Q • RSP (51 
TShIFT  ■ RSPIGI 
NT  ■ N 
C 

OATA  T-OPI/6. 25318531/ 

« ■ TWOPI #F 
IF(O.EO.OoO)  GO  TO  I 

ALPHA  » 0. 5*9/0* 

1 continue 

IFlA.NE.O.O)  30  TO  2 

TP  ■ AT  AN (•/ ALPHA) /V 

A ■ 1 .0/(EAP(-ALPHA*TP) •SIN(B«TPI 1 

2 CONTINUE 

00  10  I ■ 1 *NT 

A ■ Till  • IShIFT 

FTdl  ■ A*£X>(-ALPHA**1*SIN(B**> 

10  CONTINUE 

A • - TSHIFT 

1 • 0 

IS  I ■ I ♦ 1 

IF(Td).GE.A)  GO  TO  20 
FTdl  * 0.0 
GO  TO  15 
20  CONTtNUE 
RETURN 
ENO 

A 7/8/9  (EOR)  EXEMPT  ECL 


DEMONSTRATES  USER  SUBROUTINE  CAPABILITIES 


file  OMPSIN 

s user  defined  subroutine  to  calculate  a damped  sin  by 

S F (Tl «R2«EXP(-R3*T) •SIN(2*PI*R4#T)*U(T) 

TIME  ■ L1N128*0.5£m6  S INITIALIZE  GLOBAL  TIME  ARRAY 

S INITIALIZE  GLOBAL  REAL  ARRAY  FOR  SUBROUTINE  PUlSEl 

R2  ■ 0.  S FOR  R2*0.  THEN  F(T)  CALCULATED  SUCH  THAT  MAX  F(T)«1.0 
R3*  9.366E6  /✓  p*  ■ 1.953E7 

R5*0.  S SPECIFICS  A 0 rather  Than  a Damping  FACTOR  FOR  R2  NOT  EOUAL  TO  0. 
RA«0.  % OPTIONAL)  TINE  SHIFT  FOR  FUNCTION  SHIFT  FROM  0 

Y • PULSE  1 (TIME)  S USER  FUNCTION  PUlSEI  CALLED 

TITLE  Y • OAMPED  $IN  FUNCTION  CALCULATED  1-12-77  FROM  EXEMPT  UITm  PULSEI  • 
LIST  T 
PRTPLT  Y 
PLOT  Y 
ENO 

• G/T/8/9  (EOF)  ENO  OF  JOB 


Figure  16 


Job  Stream  and  ECL  for  Example  4 


UST  or  OAT*  SET  » 

DAMPED  SIN  rUNCTION  CALCULATED  1*12-77  FWO*  EXEMPT  KITH  PULSE1 


NOEX 

time 

VALUE  NDEX 

TINE 

V.'LUE  NOEX 

THE 

VALUE 

1 

0. 

0. 

44 

1.693E-07 

.216 

8T 

3.396E-07 

-3.063E-02 

2 

3.937E-09 

.303 

45 

1.732E-87 

.149 

88 

3.4256-07 

-4.220E-02 

3 

7.8746-09 

.959 

46 

1 . 772E-87 

5.3146-02 

89 

3. 4566-07 

-4*3586-02 

A 

1.1816-08 

.999 

47 

1 .8116-07 

-4.740e-62 

90 

3. 5346-07 

-3.519C-02 

9 

1.5756-08 

.907 

48 

1.8506-07 

-.130 

91 

3.5436-07 

-1.9596-02 

A 

1.9696-08 

.521 

49 

1. 8906-07 

-.178 

92 

3.5936-07 

-7.4706-04 

7 

2.3626-08 

.217 

50 

1.9296-07 

-.163 

93 

3.6226-07 

1.6926-02 

a 

2.7566-08 

-.207 

51 

1.9696-07 

-.  147 

94 

3.6516-07 

2.9576-02 

9 

3.1506-08 

•*  554 

52 

2. 009E-07 

-8. 128E-02 

95 

3.701E-07 

3.4776-02 

10 

3.5436-08 

-.754 

53 

2.0476-07 

-1 • 79SE-03 

96 

3.7406-07 

3.1866-02 

11 

3.937E-08 

-.772 

54 

2.0876-07 

7.2436-02 

97 

3.7906-07 

2.2126-02 

12 

4.331E-08 

-.518 

55 

2.1256-07 

.125 

98 

3.6196-07 

8.1416-03 

13 

4.7246-08 

-.337 

56 

2.1656-07 

.147 

99 

3.8586-07 

-6.6466-03 

14 

5.1186-08 

-1.394E-03 

57 

2.205E-07 

.134 

100 

3.8986-07 

-1.8916-02 

IS 

5.5126-08 

.310 

58 

2.2446-07 

9.2346-02 

101 

3.9376-07 

-2.6106-02 

16 

5.9066-08 

.531 

59 

2.2836-07 

3.3256-02 

102 

3.9766-07 

-2.6996-02 

17 

6.2996-08 

.518 

60 

2.3236-07 

-2.901E-02 

103 

4.016E-07 

-2.1836-02 

IS 

6.6936-08 

.562 

61 

2.3626-07 

-8. 041 E-02 

104 

4.0556-07 

-1.2196-02 

19 

7.0876-08 

.385 

62 

2.4026-07 

-.110 

105 

4.0946-07 

-5.2866-04 

20 

7.4806-08 

.136 

63 

2.4416-07 

-.114 

106 

4.1346-07 

1.0426-02 

21 

7 .8746-08 

-.126 

64 

2.4806-07 

-9.1456-02 

107 

4.1736-07 

1.8286-02 

22 

8.2686-08 

-.342 

65 

2.5206-07 

-5.0586-02 

108 

4.2136-07 

2.1526-02 

23 

8.6616-08 

"•  466 

66 

2.5596-07 

-1.3926-03 

109 

4.2526-07 

1.9766-02 

24 

9.0556-08 

-.478 

67 

2.5996-07 

4.4616-02 

110 

4.291E-07 

1.3736-02 

25 

9.4496-08 

-.383 

68 

2.6396-07 

7.7446-02 

111 

4.3316-07 

5.0926-03 

26 

9.8436-08 

-.210 

69 

2.6776-07 

9.0756-02 

112 

4.3706-07 

-4.0656-03 

27 

1.0246-07 

-2. 3456-03 

70 

2.7176-07 

8.2986-02 

113 

4.4096-07 

-1.1676-02 

28 

1.0636-07 

.191 

71 

2.7556-07 

5.7346-02 

114 

4.4496-07 

-1.6146-02 

29 

1.1026-07 

.328 

72 

2.7956-07 

2.0906-02 

US 

4.4986-07 

-1.6726-02 

30 

1.1426-07 

.383 

73 

2.8356-07 

-1.7756-02 

116 

4.5286-07 

-1.3546-02 

31 

1.1816-07 

.349 

74 

2.8746-07 

-4.9636-02 

117 

4.5576-07 

-7.5846-03 

32 

1.2206-07 

.239 

75 

2.9136-07 

-6.8236-02 

118 

4.6086-07 

-3.6826-04 

33 

1.260E-07 

8. 4916-02 

76 

2.9536-07 

-7.0366-02 

119 

4.6466-07 

6.4166-03 

34 

1.2996-07 

-7. 744E-02 

77 

2.9926-07 

-5.6736-02 

120 

4.6956-07 

1.1296-02 

35 

1.3396-07 

-.211 

78 

3.0316-07 

-3.1486-02 

121 

4.72»E-07 

1.3326-02 

36 

1.378F-07 

-.288 

79 

3.0716-07 

-1.034E-03 

122 

4.7546-07 

1.2246-02 

37 

l. 4176-07 

-.296 

80 

3.1106-07 

2.747E-02 

123 

4.8036-07 

8.5286-03 

38 

1.4576-07 

-.238 

81 

3.1506-07 

4.7966-02 

124 

4.8436-07 

3.1856-03 

39 

1.4966-07 

-.131 

82 

3.1896-07 

5.6176-02 

125 

4,8926-07 

-2.4866-03 

40 

1.5356-07 

-2.1 78E-03 

83 

3.2296-07 

5.1436-02 

126 

4.9216-07 

-7.2026-03 

41 

1.575E-07 

.118 

84 

3.2696-07 

3.5616-02 

127 

4.9516-07 

-9.9846-03 

42 

1.6146-07 

.203 

85 

3.3076-07 

1.3026-02 

128 

5.0006-07 

-1.0356-02 

43 

1.6546-07 

.237 

86 

3.3456-07 

-1.0966-02 

• #••4 

ENO  Of  LIST  •***•**•1 

Figure  17.  LIST  Output  of  Example  4 


98 


1.1 


J.»  *A 


A * 


.*• 


.20  * 0 


0 0 


00 


0 0 


0 000 

0 0 0 00 

0 0 00  0 00 
0 0 0 00  0000  0000  0000 
» * 0 0 0 0 00  00  00  0000  00000  0000 
0 00  00 


0 000 


-.20 


0 0 0 
* 


A A 


A A 

A 

00 


••  S.OO'rOOl.OCE-OTl.SOE-OTZ.OaE-OTZ.StE-OTI.OOE-OTl.JOE-OTA.aOE-OTA.m-OTS.OOE-OT 

TK 

sat  a ser  r 

TITLC  ■ OAMPCO  SI*  'UNCTION  CALCULATES  1-12-TT  FRSH  EXEMPT  g|TH  PULSE  l 
A - MAO  T 


Figure  18.  PRTPLT  Output  of  Example  4 
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